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Abstract. We study the susceptibility, i.e., the mean size of the com- 
ponent containing a random vertex, in a general model of inhomogeneous 
random graphs. This is one of the fundamental quantities associated to 
(percolation) phase transitions; in practice one of its main uses is that 
it often gives a way of determining the critical point by solving certain 
linear equations. Here we relate the susceptibility of suitable random 
graphs to a quantity associated to the corresponding branching process, 
and study both quantities in various natural examples. 



1. Introduction 

The susceptibility x{G) of a (deterministic or random) graph G is defined 
as the mean size of the component containing a random vertex: 

x{G) = \G\-^ (1.1) 

where C{v) denotes the component of G containing the vertex v. Thus, if 
G has n = \G\ vertices and components Cj = Ci{G), i = 1, . . . where 
K = K(G) is the number of components, then 

X{G)■.= j2—\C^\ = -Y.\C^\'. (1.2) 

i=l i=l 

Later we shah order the components, assuming as usual that |Ci| > IC2I > 

When the graph G is itself random, in some contexts (such as percolation) 
it is usual to take the expectation over G as well as over v. Here we do not 
do so: when G is random, x{G) is a random variable. 

Remark 1.1. The term susceptibility comes from physics. (We therefore 
use the notation x^ which is standard in physics, although it usually means 
something else in graph theory.) The connection with the graph version is 
through (e.g.) the Ising model for magnetism and the corresponding random- 
cluster model, which is a random graph where the susceptibility ()1.2p . or 
rather its expectation, corresponds to the magnetic susceptibility. 
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The susceptibility has been much studied for certain models in mathe- 
matical physics. Similarly, in percolation theory, which deals with certain 
random infinite graphs, the corresponding quantity is the (mean) size of the 
open cluster containing a given vertex, and this has been extensively stud- 
ied; see e.g. Bollobas and Riordan [7]. In contrast, not much rigorous work 
has been done for finite random graphs. Some results for the Erdos-Renyi 
random graphs G{n,p) and G{n,m) can be regarded as folk theorems that 
have been known to experts for a long time. Durrett proves that the 
expectation Ex{G{n,p)) = (1 — A)~^ -|- 0{l/n) if p = X/n with A < 1 fixed. 
The susceptibility of G{n,p) and G{n,m) is studied in detail by Janson and 
Luczak [22]. For other graphs, one rigorous treatment is by Spencer and 
Wormald [31], who study a class of random graph processes (including the 
Erdos-Renyi graph process) and use the susceptibility to study the phase 
transition in them. 

The purpose of the present paper is to study x{G^{i^j i'^)) for the inhomo- 
geneous random graph G^{n, k) introduced in Bollobas, Janson and Riordan 
[J]; this is a rather general model that includes G{n,p) as a special case. In 
fact, much of the time we shall consider the more general setting of 0]. We 
review the fundamental definitions from in Section [2] below. 

We consider asymptotics as n — > 00 and use standard notation such as Op, 
see e.g. 0]. All unspecified limits are as n — > 00. 

Remark 1.2. We obtain results for G{n,p) as corollaries to our general 
results, but note that these results are not (and cannot be, because of the 
generality of the model G^{n, k)) as precise as the results obtained by Janson 
and Luczak The proofs in the two papers are quite different; the proofs 
in [l^] are based on studying the evolution of the susceptibility for the 
random graph process obtained by adding random edges one by one, using 
methods from stochastic process theory, while the present paper is based on 
the standard branching process approximation of the neighbourhood of a 
given vertex. It seems likely that this method too can be used to give more 
precise results in the special case of G{n,p), but we have not attempted 
that. (Durrett [13] uses this method for the expectation Ex(G(n,p)).) 

The definition (II. 2p is mainly interesting in the subcritical case, when all 
components are rather small. In the supercritical case, there is typically one 
giant component that is so large that it dominates the sum in ()1.2p . and 
thus x{G) ~ ]Ci]^/n. In fact, in the supercritical case of [1, Theorem 3.1], 
\Ci\ = Qpin) and \C2\ = Op(n), and thus 

K K 

\Ci\^ = \Ci\^ + 0{\C2\ J2 = + 0{\C2\n) = (1 + Op(l))lCil2. 

i=l i=2 

(See also [22, Appendix A] for G{n,p).) In this case, it makes sense to 
exclude the largest component from the definition; this is in analogy with 
percolation theory, where one studies the mean size of the open cluster 
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containing, say, vertex 0, given that this cluster is finite. We thus define the 
modified susceptibility x{G) of a finite graph G by 

x(G):=^X;iC^l'- (1-3) 

i=2 

Note that we divide by n rather than by n — \Ci\, which would also make 
sense. 

In the uniform case, one interpretation of xiG) is that it gives the rate of 
growth of the giant component above the critical point. More generally, if 
we add a single new edge chosen uniformly at random to a graph G, then the 
probability that Ci becomes joined to Ci is asymptotically 2|Cj| |Ci|/n^, and 
when this happens |Ci| increases by |Cj|. Thus (under suitable assumptions), 
the expected increase in \Ci\ is asymptotically 2\Ci\ \Ci\'^/n = 2|Ci|x(G'). 

The results in ^ on components of G^{n, k) are based on approximation 
by a branching process X^, see Section [2j We define (at least when /^(5) = 1, 
see Section [2j) 

x(aj) :=E|X,| G [0,oo], (1.4) 
x(k) :=E(|X,|;|X,| < oo) G [0,cx)]. (1.5) 

Thus, x('^) = xi'^) when the survival probability p{k) := P(|3£k| = oo) = 
(the subcritical or critical case), while x('^) = oo > xi''^) when p{k) > (the 
supercritical case). 

Our main result is that under some extra conditions, the [modified] sus- 
ceptibility of G^{n, k) converges to x('^) [x('^)]i see Section [Hand in partic- 
ular Theorems 14.71 and 14. 8i 

We also study the behaviour of x(^^) and x('^^) ^ functions of the 
parameter A G (0,oo), and in particular the behaviour at the threshold for 
existence of a giant component, see Section [5l this provides a way to use 
the susceptibility to find the threshold for the random graphs treated here. 
(See, e.g., Durrett [17] and Spencer and Wormald [^] for earlier uses of this 
method.) 

Finally, we consider some explicit examples and counterexamples in Sec- 
tion [H 

Remark 1.3. We believe that similar results hold for the 'higher order 
susceptibilities' 

' ' v<^V(G) ' ' i 

but we have not pursued this. (For G{n,p), see [2^].) 

Acknowledgements. Part of this work was carried out during the pro- 
gramme "Combinatorics and Statistical Mechanics" at the Isaac Newton 
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We review the fundamental definitions from [J; but refer to those 
papers for details, as well as for references to previous work. In terms of 
motivation and applications, our main interest is the model k) of 0], 

but for the proofs we sometimes need (or can handle) different generality. 

2.1. The random graph models. In all variations we start with a measure 
space (5, /x) with < /u(5) < oo (usually, but not always, is a probability 
measure, i.e., /u(5) = 1), and a kernel on it, i.e., a symmetric non-negative 
measurable function k : S x S ^ [0, oo). We assume throughout that k is 
integrable: y) diJ.{x) d/i(y) < oo. 

2.1.1. The general inhomogenous model. To define G^(n, «;), we assume that 
we are given, for each n > 1 (or perhaps for n in another suitable index set 
Z C (0, oo)), a random or deterministic finite sequence x„ = (xi, X2, • • • , x^^) 
of points in S. (For simplicity we write Xi instead of xf^^) We denote the 
triple {Sjfi, (x„)„>i) by V and define the random graph G„ = G^{n,K) by 
first sampling x„ = {xi,X2, . . . ,Xv„) and then, given x„, taking the graph 
with vertex set {1, . . . ,Vn} and random edges, with edge ij present with 
probability min(K(xi, Xj)/n, 1), independently of all other edges. (Alterna- 
tively, and almost equivalently, see [J] and we may use the probability 
1 — exp{—K{xi,Xj)/n).) We interpret Xi as the type of vertex i, and call 
(5, /J.) the type space. 

We need some technical conditions. In 0], we assume that 5 is a separable 
metric space and fi a Borel measure; we further assume that if i^n is the 
(random) measure "^Jl^ Sx^, then f„ — ^ /u (with weak convergence 
of measures); in this case V is called a generalized vertex space. In the 
standard special case when Vn = ri and ^i{S) = 1, V is called a vertex 
space. Furthermore, in [3] it is assumed that the kernel k is graphical on V, 
which means that k, is integrable and a.e. continuous, and that the expected 
number of edges is as expected, i.e., that Ee(G^(n, K))/n — > ^ k. 

Many of the results in [4;] extend to sequences G^(n, k„), where (k^) is 
a sequence of kernels on V that is graphical on V with limit k; see [J] for 
the definition and note that this includes the case when all k„ = k for some 
graphical kernel n. 



As shown in [4 Section 8.1], if V is a generalized vertex space, we may 
condition on (x„,)„>i, and may thus assume that the x„ and, in particular, 
Vn are deterministic. Replacing the index n by Vn, and renormalizing ap- 
propriately (see Remark 12.11 below), we may reduce to the case of a vertex 
space. 



2. Preliminaries 
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2.1.2. The i.i.d. case. Another, often simpler, case of the general model is 
when is an arbitrary probability space and {xi, . . . ,Xn) are n i.i.d. 
points with distribution ;u; in this case k can be any integrable kernel. 
This case was unfortunately not treated in jl], but corresponding results 
are shown for this case (and in greater generality) in 0. In this case we 
call V = (5, /i, (x„)„>i) an i.i.d. vertex space. In this case, to unify the 
notation, a graphical kernel is thus any integrable kernel. Many results for 
this case extend to suitable sequences of kernels, for example assuming that 
ll'^n — — > 0, as then the general setting below applies. 

2.1.3. Cut-convergent sequences. To define the final variant we shall con- 
sider, we briefly recall some definitions. (A variant of) the Frieze-Kannan [l~ 
cut norm of an integrable function : 5^ ^ M is simply 



sup / f{x)W{x,y)g{y)dfi{x)dfi{y). 

ll/lloc,l|<?||oc<1^52 

Given an integrable kernel k and a measure-preserving bijection r : 5 — > 5, 
let K^"^^ be the corresponding rearrangement of n, defined by 

K^''\x,y) = K{T{x),T{y)). 

We write k ~ k' if k' is a rearrangement of k. Given two kernels k, k' on 



[0,1], the cut metric of Borgs, Chayes, Lovasz, Sos and Vesztergombi [lO(] 
may be defined by 

^□(k,k')= inf — (2-1) 
There is also an alternative definition via couplings, which also applies to 



kernels defined on two different probability spaces; see |ld: Iq. 

Suppose that A„ = (aij) is an n-by-n symmetric matrix with non-negative 
entries; from now on any matrix denoted An is assumed to be of this form. 
Then there is a random graph G„ = G{An) naturally associated to A^. the 
vertex set is {1,2, . . . ,n}, edges are present independently, and the proba- 
bility that ij is an edge is min{ ajj/n, 1}. Given An, there is a corresponding 
kernel ka„ on [0, 1] with Lebesgue measure: divide [0, 1]^ into squares of 
side 1/n in the obvious way, and take the value of ka„ on the (i, j)th square 
to be Qij. Identifying An and the corresponding kernel, as shown in [q], 
many of the results of [3] apply to G„ = G{An) whenever 6\j{An, k) — > for 
some kernel k on [0, 1] (or, more generally, on some probability space S). 

If An is itself random, then G{An) is defined to have the conditional dis- 
tribution just described, given A„. Any results stating that if S\j{An, k) — > 
then G{An) has some property with probability tending to 1 apply also if 
(An) is random with 5\j{An, k) — ^ 0. (One way to see this is to note that 
there is a coupling of the distributions of the An in which 6^{An, n) ^ 
a.s., and we may then condition on (An).) 

Moreover, as shown in ^, Sections 1.2 and 1.3], such results apply to 
the models described in the previous subsections, since in each case the 
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(random) matrices of edge probabilities obtained after conditioning on the 
vertex types converge in probability to k in 6]j. 

2.2. The corresponding branching process. Given an integrable kernel 
K on a measure space (5,/u), let X^ix), x £ S, he the multi-type Galton- 
Watson branching process defined as follows. We start with a single particle 
of type X in generation 0. A particle in generation t of type y gives rise to 
children in generation t + 1 whose types form a Poisson process on S with 
intensity K{y , z) diJ,{z) . The children of different particles are independent 
(given the types of their parents). 

If is a probability measure, we also consider the branching process 
defined as above but starting with a single particle whose type has the 
distribution ^u. 

Let |Xk(x)| denote the total population of Xk{x), and let 

Pk{K;x) ■.= F{\X^{x)\ = k), A; = l,2,...,oo, (2.2) 

and 

Pk{n) := j pk{K;x)dfj,{x), k = 1,2, . . . ,oo. (2.3) 

Thus, when /u(5) = 1, Pki^) is the probability P(|3Ck| = k). 
For convenience we assume that 

j^K{x,y) dn{y) <oo (2.4) 

for all X G S; this implies that all sets of children are finite a.s. This is 
no real restriction, since our assumption that Jg2 k < oo implies that (j2.4p 
holds for a.e. x, and we may impose (j2.4p by changing «; on a null set, which 
will a.s. not affect X^- (Alternatively, we could work without ()2.4p . adding 
the qualifier "for a.e. x" at some places below.) 

Since a.s. all generations of X^i^x) are finite, it follows that poo{n;x), the 
probability that the branching process is infinite, equals the survival proba- 
bility of Xk{x), i.e., the probability that all generations are non-empty. We 
use the notation p{k; x) := Poo{k; x); for typographical reasons we sometimes 
also write Pk.{x) = p{n;x). Similarly, we write p{n) := Poo{k); if p{S) = 1, 
this is the survival probability of X^- 

We are interested in the analogue of the mean cluster size for the branch- 
ing processes. For Xk{x), we define 

x(k;x) :=E(|X,(x)|) = Yl kpk{K;x), (2.5) 

l<fe<oo 

x{k;x) ■.= E[\Xf,{x)\;\Xf,{x)\ < oo) = ^ kpk{K;x); (2.6) 

l<fc<oo 
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thus x{f^'ix) = xi^^'ix) < OO if p{k;x) = 0, and x{f^]^) ^ xi^^'j^) = oo if 
p{k; x) > 0. Further, let 



x{K):=fi{S)-' x{^;x)dfi{x)=fi{S)-' ^ kpk{K), (2.7) 

•^^ l<k<oo 

xiK):=piS)-' j^xi'^;x)dfi{x)=piS)-' ^ kphi^). (2.8) 



l<fc<oo 



Thus, if p{S) 



X(M=E(|X,|), (2.9) 
X(M =E(|X,|;|XJ < oo). (2.10) 

Remark 2.1. For a generaUzed vertex space, where p{S) may differ from 
1, we may renormalize by replacing p and k by 

p := p{S)^^p and k' := p{S)k. (2-11) 

This will not affect Xk,{x), and thus not x('«; 2;) and x(«^; x); further, because 
of our choice of normalization in ()2.7p and (12. Sp . x('^^) and x('^^) also remain 
unchanged. Hence, results for generalized vertex spaces follow from the case 
when p{S) = 1. 

2.3. Integral operators. Given a kernel k on a measure space {S,p), let 
be the integral operator on (5, p) with kernel k, defined by 

{TJ){x):= [ K{x,y)f{y)dp{y), (2.12) 

for any (measurable) function / such that this integral is defined (finite 
or +00) for a.e. x. (As usual, we shall assume without comment that all 
functions considered are measurable.) Note that T^/ is defined for every 
/ > 0, with < T^f < 00. 
We define 

IIT^II := sup{||T,/||2 : / > 0, ||/||2 < l} < 00. (2.13) 

When finite, \\Tf^\\ is the norm of as an operator in L^{S, p). We denote 
the inner product in (real) L'^{p) by {f,g) = {f,g)^ '■= J^fgdp, and the 

norm by II/II2 := (/, 

One of the results of [J] is that the function p^ix) = p{k; x) is the unique 
maximal solution to the non-linear functional equation 

/ = 1 - e-^-f, / > 0. (2.14) 

Moreover, if \\Tf^\\ < 1, then p^ = and thus p{k) = 0, while if ||Tk|| > 1, 
then > on a set of positive measure and thus p{k) > 0. (This extends 
to generalized vertex spaces by the renormalization in Remark 12. Ij note 
that 7k and \\Ti^\\ are not changed by the renormalization.) The three 
cases IITkII < 1, \\Tf^\\ = 1 and ||Tk|| > 1, are called subcritical, critical and 
supercritical, respectively. 
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Given a kernel «; on a type space (5, let /2 be the measure on S defined 

by 

d/2(x) := (1 -p(K;x))d^(x). (2.15) 
(This is interesting mainly when k is supercritical, since otherwise fl = fi.) 
The dual kernel k is the kernel on (5, /I) that is equal to function. We 

regard as an operator acting on the corresponding space L'^{fi). Then 
||r^|J< 1; typically \\T-^\\ < 1 when k is supercritical, but equality is possible, 
see [J, Theorem 6.7 and Example 12.4]. 
Note the explicit formula 



{TTif){x) := / K{x,y)f{y)djl{y) = / K{x,y)f{y){l-p{K;y))dfi{y), 
Js Js 

(2.16) 

i.e., T^f = T«;(/(l - Pk))- Note also that 

JUS) = ^(1 - ^)) d/u(x) = p(5) - pin); (2.17) 
if /i(5) = 1, this is the extinction probability of X^- 

2.4. Small components. Let Nk{G) denote the number of vertices in com- 
ponents of order /c in a graph G. (Thus the number of such components is 
Nk{G)/k.) We can write the definition (II. 2|) as 

x(G) = pE^^ =E^^- (2.18) 

By |3, Theorem 9.1], if (k n) is a graphical sequence of kernels on a vertex 
space V with limit k and G„ = G^{n, Kn), then, for every fixed A; > 1, with 
N>k := T,j>k^j and /9>fc := Efc<i<ooPi> we have 

iV>fc(G'„)/n^p>fe(K), (2.19) 

and thus 

Nk{Gn)/n^Pk{K). (2.20) 

This extends to generalized vertex spaces by normalization (if necessary first 
conditioning on (x„)„>i) as discussed in y, Subsection 8.1]. Furthermore, 
(|2.2U|) holds also on an i.i.d. vertex space for a constant sequence k„ = k, 
with K integrable, by [1, Lemma 21]. 

Even more generally, by 0, Lemma 2.11], the same conclusions hold when 
Gn = G{An) with 5g(j4„,k) — > 0, and hence when G„ = G{An) with 
5\j{An, k) — ^ 0; this implies the two special cases above. 

2.5. The giant component. If k is irreducible (see Q] for the definition), 
then under any of our assumptions we have 

\Ci{Gn)\/n^ p{k) (2.21) 

and 

\C2{Gn)\/n^0; (2.22) 



SUSCEPTIBILITY IN INHOMOGENEOUS RANDOM GRAPHS 9 

see [1, Theorems 3.1 and 3.6] or 0, Theorem 1.1]. 

2.6. Monotonicity. We note a simple monotonicity for x; there is no cor- 
responding result for x- 

Lemma 2.2. If H is a subgraph of G with the same vertex set, then x{H) < 
X{G). 

Proof. Immediate from the definition (jl.l[) . □ 

3. Branching processes 

For branching processes, as is well-known, the mean cluster size can be 
expressed in terms of the operators and T-fi. We write 1 for the constant 
function 1 on 5. 

Lemma 3.1. For any integrable kernel k on a type space we have 

oo 

x{K;x)=Y,ni{x), (3.1) 
j=0 

oo „ oo 

x{k) = f,{S)-' / Til{x) d/i(x) = Y.^Til, 1)^, (3.2) 

j=Q j=0 

oo 

x{K;x) = {l-p{K;x))YTil{x), (3.3) 

j=Q 

oo „ oo 

/i(5)-i J Tllix) dKx) = KSr' Yi^ih (3.4) 



j=o -^-s j=0 



Proof. Let fj{x) be the expected size of generation j in Xk(x). Then, for 
every j > 0, by conditioning on the first generation, 

fj+i{x) = fj{y)n{x,y) dfi{y) = T^fjix), 

and thus, by induction, fj = T^fo = T^l. Hence, ()3.ip follows by summing. 
Recalling the definition (j2.7p . relation (j3.2p follows immediately. 

It is easy to see that if we condition Xi^{x) on extinction, we obtain another 
similar branching process Si^ix) with /i replaced by p.. Hence, is replaced 
by T^, and (j3.3p follows from 

E{\X^{x)\;\X^{x)\ < oo) = (1-p(k;x))E(|X^(x)| | |X^(x)| < oo) 

= (1-p(k;x))E(|X,(x)|) 
and (f3T]) . Finally, (jS^i]) follows by ([2:8|) and integration, recalhng ([2J5]) . □ 
Often, it is convenient to assume for simplicity that /u(5) = I. 



10 



SVANTE JANSON AND OLIVER RIORDAN 



Lemma 3.2. Let n he an integrable kernel on a type space {S, fi) with fJ-{S) = 
1 . Then 

oo 

x{k) = Y^iTll, 1)^ = KS)X{^) = (1 - Pii^m^- 

j=0 

Proof. Use (I3.4p for k and fi and (13. 2p for k and /I, together with (I2.17p . □ 

Theorem 3.3. Let k be an integrable kernel on a type space (5,^) with 
fi{S) = 1. 

(i) If K is subcritical, i.e., \\Ti^\\ < 1, thenx{i^',x) = (/ — Tfj)"-*^! a.e., and 

= ((/-r,)-H,i)^ < oo. 

(ii) Suppose that k is supercritical, i.e., \\Ti^\\ > 1, and also that \\Tfi\\ < 1. 
Thenx{K;x) = {I- p,){I -T^^yH a.e., andxi^) = {{I-Tf.)-H,l)ji < 
oo. 



The conditions of (ii) hold whenever \\Tf^\\ > 1, k is irreducible, and J^2 



OO. 



Proof. An immediate consequence of Lemma 13.11 since in these cases the 
sumsX]jlo^« = (^-^k)'-^ andX^°^or| = {I-Tf=)~^, respectively, converge 
as operators on L'^[p) and L'^{'j2). For the final statement we use [1, Theorem 
6.7], which yields \\Tfi\\ < 1. □ 

In fact, for the last part one can replace the assumption that J^2 < oo 
by the weaker assumption that is compact; this is all that is used in the 
proof of [3, Theorem 6.7]. 

In the critical case, when \\Tfi\\ = 1, we have x('*) = x(^)- We typi- 
cally expect the common value to be infinite, but there are exceptions; see 
Section I 



Theorem 3.4. (i) If k is critical and is a compact operator on L'^{ii), 
then = oo. In particular, this applies if i^{x.,y)'^ d/i(x) dij,{y) < oo. 
(ii) If K is supercritical, then xi^^) = oo- 

Proof, (i): If J^2 < oo, then is a Hilbert-Schmidt operator and thus 
compact. 

is always self-adjoint (when it is bounded), so if is compact and crit- 
ical, then it has an eigenfunction ip with eigenvalue \\Tf^\\ = 1; moreover, the 
eigenspace has finite dimension and there is at least one such eigenfunction 
tpi > (with 11^1 II 2 = 1, say), see Lemma 5.15 in Q] and its proof, where 
only compactness is used. There may also be eigenfunctions with eigenvalue 
— 1, so we consider the positive compact operator T^ and let V'l, • • • > V'm be 
an orthonormal basis of the eigenspace for the eigenvalue 1 of T^. The or- 
thogonal complement is also invariant, and T^ acts there with norm R < 1. 
Hence, 

m m 

(T^l, 1) = ^(1, + 0{Rn - 5^(1, ^^)'. 

i=l i=l 
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Since the terms in the sum are non-negative and (1, "^i) = J ipidfi > 0, the 

hmit is strictly positive and thus YlJLoiT^lA) cannot converge. Since the 

terms in this sum are non-negative, (j3.2p yields x(k) = n{S)~^ Yl'jLoiTK^, 1) = 
oo. 

(ii): By [4, Theorem 6.1] we have P(|3Ck| = oo) = p{k) > 0, so x(^) = 
oo. □ 

In the subcritical case, we can find xi^^) by finding (/ — T^)"^!, i.e., by 
solving the integral equation f = T^f + 1. Actually, we can do this for any 
K, and can use this test of whether xi^^) < oo. 

Theorem 3.5. Let k he a kernel on a type space Then the following 

are equivalent: 

(i) X(«;) < oo- 

(ii) There exists a function f > in L^i^n) such that (a.e.) 

f = Tf + l. (3.5) 

(iii) There exists a function f > in L^{lj) such that (a.e.) 

f>Tf + l. (3.6) 

When the above conditions hold, there is a smallest non-negative solution f 
to (13. 5p . that is also a smallest non-negative solution to (13. 6p ; this minimal 
solution f equals xii^i x), and thus x(^) = 1^{'S)~^ /d/^. 

Proof. Recalling (j3.ip . let g{x) := xi^^'^x) = Yl'jLoTKl{x); this is a function 
S — > [0, oo] with T^g = X^^^i T^l = g — 1, so g satisfies both (13. 5p and (|3.6p . 
Further, J^gd^i = ij.{S)x{k) by (|3.2p . Hence, if (i) holds, then g S L^{fJ.); 
consequently, g satisfies (ii) and (iii). (Note that then g is finite a.e.) 
Conversely, if / > solves (]3.5p or (13. 6j) . then, by induction. 



n-l 

f>Y,ni+T:f 

j=0 

for every n > 1. Thus / > E"=o^«l> and letting n ^ oo yields f > g. 
Hence, if (ii) or (iii) holds, then g £ L^(fi), and (i) holds. Further, in this 
case, f > g, which shows that g is the smallest solution in both (ii) and (iii), 
completing the proof. □ 

Note that in the subcritical case, (13. 5|) always has a solution in L^(//); 
cf. Theorem 13.31 In Section 16.31 we give an example where k is critical 
and (13. 5p has a solution that belongs to L^ifJ^), but not to L'^{p). (We do 
not know whether there can be a non- negative solution in L'^{fJ.) with k 
critical.) Moreover, in this example, both in subcritical and critical cases, 
there is more than one non-negative solution in L^(^). However, we can 
show that there is never more than one non-negative solution in L^(^). 
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Corollary 3.6. Suppose that there exists a function f > in Lp'{^i) such 
that (13.51) holds. Then f is the unique non-negative solution to (13. 5p in 



L^ifJ-), xiK;x) = f{x) and x{k) = fJ-iS) ^ f dfi. 

Proof. Let g be the smallest non-negative solution, guaranteed to exist by 
Theorem [331 and let h = f - g > 0. Since < h < f , h £ L'^iji). Then 
Th = Tf-Tg={f-l)-{g-l)= h, and 

(/, h) = {Tf + l,h) = {Tf, h) + (1, h) = if, Th) + (1, h) = if, h) + (1, h). 
Hence = = J /id^, so /i = a.e., and f = g. □ 



We begin with a general asymptotic lower bound for the susceptibility. 
This bound depends only on convergence of the number of vertices in compo- 
nents of each fixed size, so it applies under any of the assumptions described 
above. More precisely, we state the results in the setting of Subsection l2.1.3| 
as noted there they then apply (by conditioning) to G^(n, k„) under the as- 
sumptions in Subsection 12.1.11 or Subsection 12.1.21 As usual, we say that G„ 
has a certain property with high probability, or whp, if the probability that 
Gn has this property tends to 1 as n — > oo. 

Recall that a matrix denoted An is assumed to be symmetric, n-by-n and 
to have non- negative entries. 

Theorem 4.1. Let k be a kernel and {An) a sequence of (random) matri- 
ces with d\j{An,K) — ^ 0, and set Gn = G{An). Alternatively, let Gn = 
G^{n,Kn) satisfy the assumptions of Subsection [Kl.l\ or Subsection \2. 1 
Then, 

(i) for every b < x(^); whp x{Gn) > b, and 

(ii) for every b < x(^)) whp xiGn) > b. 
Moreover, liminf Ex(Gn) > xi^^) o'^^ liminf Ex(G„) > 

Proof. As noted in Subsection 12.1.31 after reducing to the vertex space case 
if necessary (and so assuming without loss of generality that fJ.{S) = 1) it 
suffices to consider the case Gn = G{An). 



(i): Let K he a fixed positive integer. Then, by (I2l8]) . (f2T9]l and (1212(1 . 



4. Main results 




k=l 




K-l 



K-l 



^ kpkjK) + Kp>k{h) = {kAK)pk{K). 
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As K ^ oo, the right-hand side tends to x('^) by monotone convergence 
and ()2.7p : hence we can choose a finite K such that the right-hand side is 
greater than b, and (i) follows. 

(ii) : By (11. 2[) and (|1.3p . if Ci is the largest component of Gn and \Ci\ > K, 
then 

X{Gn) > > K . 

n 

k=l 

On the other hand, if \Ci\ < K, then 
Hence, in both cases, using (j2.20p again, 

> f; - — ^ E fcp.(-). (4.1) 

k=l k=l 

As — > oo, the right-hand side tends to x(^)) thus we can choose K 
such that it exceeds b, and (ii) follows. 

(iii) : An immediate consequence of (i) and (ii). □ 

We continue with a simple general probability exercise. 

Lemma 4.2. Let X„ be a sequence of non-negative random variables and 
suppose that a E [0, oo] is such that 

(i) for every real b < a, whp X„ > b, and 

(ii) limsupEX„ < a. 

p 

Then Xn — > a and EX„ — > a. Furthermore, if a < oo, then Xn — > a, i.e., 
E|X„ - a| ^ 0. 

Proof. If a = oo, (i) says that X„ — ^ oo; this implies liminfEX„ > b for 
every b < oo, and thus EX„ — > oo. 

Assume now that a < oo, and let e > 0. Then, for every fe < a, by (i), 

E(X„ - a) > eF{Xn > a + e) - {a - b)F{a + e > X^ > b) - aP(X„ < b) 

> eF{Xn >a + e)-{a-b)- o(l). 

Hence 

limsupE(X„ — a) > elimsupP(X„ > a + e) — {a — b) 
and thus, since 6 < a is arbitrary, 

limsupE(X„ — a) > elimsupP(X„ > a + e). 

Since limsupE(X„ — a) < by (ii), this yields limsupP(X„ > a + e) = 
for every e > 0, which together with (i) yields Xn a. 
Moreover, the same argument yields, for every e > 0, 

liminfE(X„ - a) > eliminfP(X„ > a + e). 

Taking e = we obtain liminfEX„ > a, which together with (ii) yields 
KX„^a. □ 
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The idea is to use Lemma 14.21 with = x{Gn) and a = x('*) oi' = 
x{Gn) and a = xl^^); then condition (i) is satisfied by Theorem l4.lt and we 
only have to verify the upper bound (ii) for the expected susceptibihty. For 
convenience, we state this exphcitly. 



Lemma 4.3. Let k and Gn he as in Theorem \^.l\ 

(i) // limsupExlCn) < x{f^), then x{Gn) ^ x{>^) and ¥.x{Gn) 

(ii) // limsupExlCn) < x{i^), then x{Gn) ^ x{i^) and Ex(G„) 

xii^-)- 

Proof. By Theorem 14.11 and Lemma 14.21 as discussed above. □ 

Sometimes we can control the expectation only after conditioning on some 
(very likely) event. This still gives convergence in probablity. 



Lemma 4.4. Let k and Gn be as in Theorem \4-l\ and let £n be an event 
(depending on Gn) such that £n holds whp. 

(i) // limsupE(x(G„);f„) < x(k), then x{Gn) ^ x(«)- 

(ii) // limsupE(x(G„);f„) < x{i^), then x{Gn) x{i^). 

Proof. After conditioning on we still have Nk{Gn)/n — ^ Pk{i^) for each 
fixed /c, which is all that was needed in the proof of Theorem 14. 1[ Letting 
(/9 = X or X) since E((^(G„) | <f„) ~ E(93(G„); <5„), under the relevant as- 
sumption Lemma 14.21 tells us that the distribution of ip{Gn) conditioned on 
£n converges in probability to v'(k). But then the unconditional distribution 
converges in probability. □ 



We begin with a trivial case, which follows immediately from Lemma [4.3 
Theorem 4.5. Let k and Gn he as in Theorem\4.1\ 



(i) //x(^) = oo, then x{Gn) — ^ oo and Ex(Gn) — > oo. In particular, 
this holds if k is critical and is compact, or if k is supercritical. 

(ii) If xin) = oo, then xiGn) — ^ oo and Ex(G„) oo. 

Proof. The extra conditions in Lemma 14.31 are vacuous. For (i), we use also 
Theorem [M □ 

One way to obtain the required upper bound on the susceptibility is by 
counting paths. Let Pi = Pt{G) denote the number of paths vqVi . . .vi of 
length i in the graph G. 

Lemma 4.6. Let G be a graph with n vertices. Then xiG) < YIT^q Pe{G)/n. 

Proof. For each ordered pair {v,v') of vertices of G with v and v' in the 
same component, there is at least one path (of length > 0) starting at v and 
ending at v' . Thus, counting all such pairs, l^iP ^ Yl'^o^i- ^ 
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So far our arguments relied only on convergence of the number of vertices 
in components of a fixed size k, and so apply in very great generality. Unfor- 
tunately, bounding xiG) from above, via Lemma 14.61 or otherwise, involves 
proving bounds for all k simultaneously. These bounds do not hold in gen- 
eral; we study two special cases where they do in the next two subsections. 



4.1. Bounded kernels on general vertex spaces. In this section we 
consider G„ = G^{n,Kn), where is any uniformly bounded graphical 
sequence of kernels on a (generalized) vertex space V with limit k. In fact, 
we shall be consider the more general situation where G„ = G{An) for some 
sequence (An) of uniformly bounded (random) matrices with 6\j{An,K) — ^ 
0. From the remarks in [il, the graphs G^{n, Kn) are of this form. Note that 
this is the setting in which the component sizes were studied by Bollobas, 
Borgs, Chayes and Riordan [3]. 

Theorem 4.7. Let n he a kernel and {An) a sequence of uniformly hounded 
matrices with ^□(yl„, k) -^—^ 0, and set Gn = G{An). Alternatively, let Gn = 
G^{n,Kn) satisfy the assumptions of Subsection \2.Ll\ or Subsection \2.1.S\ 
with the Kn uniformly hounded. 

(i) We have x{Gn) ^ x{i^)- 

(ii) If K is irreducible, then xiGn) x('^)- 

The boundedness assumption is essential unless further conditions are 
imposed; see Example 16.91 The extra assumption in (ii) is needed to rule 
out the possibility that there are two or more giant components, living in 
different parts of the type space. 

Proof. As noted above, the case of a generalized vertex space V may be 
reduced to the case of a vertex space by conditioning and renormalization, 
see Subsection 12.1.11 and Remark 12.11 and the vertex space case in Subsec- 
tion 12.1.11 or Subsection 12.1.21 is a special case of the version with matrices 
An, so it suffices to consider the latter version. In particular, we may assume 
that /i(5) = 1. 

Coupling appropriately, we may and shall assume that (5n(A„, k) 0. It 
is easily seen that this and the uniform boundedness of the An imply that 
K is bounded. 

For (i), suppose ffi'st that ||Tk|| > 1. Then, since is compact, by 
Theorem 13.41 we have x('^) = c>o, and by Theorem 14. 51 we have x{Gn) oo 
as required. 

Suppose then that ||Tk|| < 1. Let k„ = ka„ denote the piecewise constant 
kernel corresponding to An- Then, letting 1 denote the vector (1,...,1), 
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and writing An = (o-"^), we have 



a). , 



io,...J<=l i=l 

= nE T] Knixi-i,Xi)dn{xo) ■ ■ ■ dn{xi) 

= n{Tij,l)^. (4.2) 

Recall that Kn and k are uniformly bounded, and 5 n(K n, n) — > 0. As noted 
in or by the Riesz-Thorin interpolation theorem [isl. Theorem VI. 10. 11] 
(for operators — > and ^ L°°), it is easy to check that this implies 
II^Knll ~^ II^kII- i^^ fact, the normalized spectra converge; see [ll|].) Since 
||Tk|| < 1, it follows that for some 5 > we have < 1 — 5 for n large 

enough, so YlA'^Kn'^^ 1)/^ — J2e II^Knll^ converges geometrically. 

For a fixed £, and kernels k, k' bounded by M, say, it is easy to check that 
|(r5l,l)^-(r^l,l);,| < IM^-'^Wk' - K\\n (see, for example, 0, Lemma 2.7]). 
Since {T^/l, 1)^ is preserved by rearrangement, we may replace — K||n by 
6\j{k' , k) in this bound. Hence, for each i, we have (T^^l,l)^ (r^l,l)^. 
Combined with the geometric decay established above, it follows that 

oo 

By Lemma 14.61 and ()4.2p we thus have 

^ oo oo 

limsupEx(G'„) < limsup - ^EP^(G„) < limsup^(T^^l, 1)^ = x('«), 

£=0 e=o 

which with Lemma l4.3r i) gives xiGn) — ^ x('^) s-s required. 

We now turn to to the proof of (ii). If ||Tk|| < 1, then p{k) = 

and x('^) = x('^)- On the other hand, x(G'n) < x(Gn), so the bound above 
gives limsupEx(G„) < xif^) = x(^); Lemma 14.3( 11) gives the result. 

Now suppose that ||Tk|| > 1. Let G„ be the graph obtained from G„ by 
deleting all vertices in the largest component Ci, and let n be the number of 
vertices of Gn- By the duality result of [13] (see also d. Theorem 12.1] for 
the case Gn = G^ {n, Kn)), there is a random sequence (Bn) of matrices (of 
random size h x h) with 5\j{Bn,K) — ^ 0, such that Gn may be coupled to 
agree whp with G{Bn)', here H := k' is k renormalized as in (j2.1ip . (Recall 
that K is regarded as a kernel on (5,/2), where /i defined by (|2.15p is not a 
probability measure.) By Remark 12.11 x(^) = xi^)- 

Note that 

\Gn\ IT' — \Cl\ P / \ / . ON 

= > 1 - p{k) (4.3) 

n n 
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by ()2.2ip . After conditioning on the number of vertices of Gn and the 
matrices Bn, we can apply part (i) to conclude that 

xiGn) = xiG{Bn)) + Op(l) ^ x(^) = Xi^)- (4.4) 
Finally, if {Cj}j>i are the components of G„, then {Cj}j>2 are the compo- 
nents of Gn, and thus by ((LS]), ((L2l), dH]), (g^]) and Lemma [32] 

X G„ = = y [1- p{k))x{k) = x{i^)- □ 

n n 

4.2. The i.i.d. case. 

Theorem 4.8. Let k be an integrahle kernel on an i.i.d. vertex space V. 
Then x{G^ {n, k)) xii^) and Ex{G^ {n, k)) x{i^)- 

Proof. Similarly to the estimate in the proof of Theorem 14.71 for any £, the 
expected number E Pi of paths of length i is 

n---{n-£) / TTminf ' , 1 ) d/i(xo) • • • d^(x^) 

Js^+^fJ[ \ n ) 

I 

<n TTK(xi_i,Xi)d^(xo) • • • d//(2;£) = n(r^l, 1)^. 

Summing over all i > 0, we see by ()3.2p that the expected total number of 
paths is at most nxin). Hence, by Lemma 



Ex{G^{n,K))<E^Pe/n<x{K). (4.5) 

The result follows by Lemma 14.31 □ 

Our next aim is to prove a similar result for x- Unfortunately, we need 
an extra assumption. We shall assume that is compact, though any 
condition guaranteeing (|4.23p below will do. 

Theorem 4.9. Let k be an irreducible, integrable kernel on an i.i.d. vertex 
space V with ||Tk|| > 1, and let Gn = G^{n,K). If is compact, then 

We do not know whether compactness, or some similar assumption, is 
necessary for this result. 

The main idea of the proof is to count the expected number of paths P 
such that P is not joined to a large component of G„ — P. We start with 
a few preparatory lemmas that hold under more general conditions than 
Theorem 14.91 itself. 

Recall that Ci = Ci(Gn) CI [n] denotes the (vertex set of) the largest com- 
ponent of Gn- As in [j], given G„, let denote the empirical distribution 
of the types of the vertices in Ci(G„), so for j4 C 5 we have 

ul,{A) = n-'\{i e CiiGn) : Xi e A}\. 
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Lemma 4.10. Let k be an irreducible, integrable kernel on an i.i.d. vertex 
space V = {S,fi, (x„,)n>i), and let A be a measurable subset of S. Then 

t^li^) IJ-K^A) := / p{K;x)dfi{x). 
J A 

More precisely, the convergence is uniform in A: given any e > there is 
an riQ such that for all n > uq and all measurable A we have 

F{\ul{A)-f^^{A)\>e)<e. 

Note that the first statement corresponds to Theorem 9.10 of [4], but, due 
to the different conditions, is not imphed by it. 

Proof. It suffices to prove the second statement. Fix e > once and for all, 
and choose ko so that p>fco('^) < +e/6; this is possible since p>k{i^) \ 
p{k) as /c — > cxD. 

We start by considering components of a fixed size. Let Nk{A) denote 
the number of vertices i of Gn such that i is in a component of order k 
and Xi ^ A. If k is bounded, then using the local coupling argument in [1, 
Section 3] it is easy to check that for each k we have Ni^{A)/n — ^ Pfc(^) := 
J^Pk{x) dfi{x), uniformly in A. Using the fact that adding or deleting an 
edge from a graph G changes the set of vertices in components of size k in 
at most 2k places, and arguing as in [3], the same statement for general k 
follows easily. 

Summing over k < ko, we thus have N<f:g{A)/n — ^ p<^ko{^)- In partic- 
ular, 

P(|iV<fc„(^)/n-/.<,„(^)| >e/5) <e/3 (4.6) 

for all large enough n and all measurable A. 

By a medium component of G„ we mean any component of size greater 
than ko other than Ci(G„). Let M denote the number of vertices in medium 
components, and M{A) the number with types in A. Since Nf:{Gn)/n — ^ Pk 
for each k and |Ci(G„)|/n — ^ p('^), we have M{Gn)/n — ^ p>kQ+i{i^) — 
p{k) < e/6. Hence, whp 

supM(A) = M{Gn) < en/5. (4.7) 

A 

Let ij={A) denote the number of vertices with types in A. Then #(vl) has 
a binomial distribution with parameters n and p{A), so for n large enough 
we have 

P(|#(^)/n-^(A)| >e/5) <£/3 (4.8) 

for all A. Finally, let Gi{A) = nv}^{A) denote the number of vertices in 
Ci{Gn) with types in A. Then 

Ci(A) = #(A) - N<kM) - M{A) + 0(1), (4.9) 

with the final 0(1) correction term accounting for the possibility that |Ci(G„)| < 
ko, so the 'giant' component is 'small'. 
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Combining equations ()4.6|) ~ ()4.9p . we see that 

¥{\CM)/n- {^l{A) - p<k,{A))\>Ae/b) <e 
for all large enough n and all A. But 

oo 

fc=fco+l 

The sum above is at least but, by choice of ko, at most e/6, so fJ.{A) — 
p<ko{A) is within e/6 of Pk{A) and the result follows. □ 

In P, Theorem 1.4], it was shown (in a slightly different setting) that 
stability of the giant component under deletion of vertices implies that the 
distribution of the size of the giant component has an exponential tail. Parts 
of this argument adapt easily to the present setting. 

First, Lemma 1.7 of [Q] shows that if k is a kernel, then the n-by-n matrices 
obtained by sampling k at i.i.d. points xi, . . . ,Xn converge in probability to 
K, with respect to the cut norm. This implies that all results of 0] asserting 
that a certain conclusion holds whp apply to the corresponding random 
graphs (see Remark 1.5]). In particular, Theorem 1.3 of [6] implies the 
following result. 

Theorem 4.11. Let k he an irreducible, integrable kernel on an i.i.d. vertex 
space V, and let Gn = G^{n, k). For every e > there is a 6 > such that 
whp we have 

pi^)-e<\CiiG'J\/n<piK) + e 

for every graph G'^ that may be obtained from Gn by deleting at most 6n 
vertices and their incident edges, and then adding or deleting at most 6n 
edges. □ 

Using this result, it is easy to get our exponential lower tail bound. Un- 
fortunately, there is a minor complication, due to the possible (but very 
unlikely) non- uniqueness of the giant component. 

Let Ci{A) = Gi{A;Gn) denote the maximum over components C of Gn 
of the number of vertices of C with types in A, so Ci{A) is within |C2(G„)| 
oiCi{A) =niy^{A). 

Lemma 4.12. Let k be an irreducible, integrable kernel on an i.i.d. vertex 
space V = (S, p, {xn)n>i) with ||Tk|| > 1, and let e > 0. Then there is a 
c = c(k, e) > such that for all large enough n, for every subset A of S we 
have 

P(C'i(A; Gn) < [pM) - e)n) < e-^^ (4.10) 

Proof. Fix A. Given a graph G on [n] where each vertex has a type in S, let 
D{G) = Da{G) be the minimum number of vertices that must be deleted 
from G so that in the resulting graph G' we have 

CM-.G')<{pM)-e)n, (4.11) 
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SO our aim is to bound P(D(G„) = 0). By Lemma 14.101 whp Ci(G„) has at 
least {fj-KiA) — e/2)n vertices with types in A. Also, by Theorem 14. Ill there 
is some 5 > such that whp deleting at most 5n vertices of G„ removes less 
than en/2 vertices from the (whp unique) giant component. It follows that 
E,D{Gn) > 6n/2 for n large; moreover, this bound is uniform in A. 

Since the condition (14. lip is preserved by deleting vertices, if G" is ob- 
tained from G by adding and deleting edges all of which are incident with one 
vertex i, and also perhaps changing the type of i, then \D{G) — D{G")\ < 1. 
We may construct G„ by taking independent variables xi, . . . ,Xn and {yij : 
1 < i < j < n} all of which are uniform on [0, 1], and joining i to j if and 
only if yij < K{xi,Xj)/n. Modifying the variables in Sj = {xj}U{yij : i < j} 
affects only edges incident with vertex j. Considering the values of all vari- 
ables in Sj as a single random variable Xj, we see that D{Gn) is a Lipschitz 
function of n independent variables, so by McDiarmid's inequality [2g] we 
have 

¥{D{Gn) = 0) < e-2(^^(G"))Vn < e-5^n/2^ 

completing the proof. □ 

It would be nice to have an exponential bound on the upper tail of the 
number of vertices in 'large' components. Unfortunately, the argument in 0] 
does not seem to go through. Indeed, the corresponding result is false in this 
setting without an additional assumption: it is easy to find a k, for which 
there is a small, but only polynomially small, chance that some vertex v has 
degree of order n. In this way one can even arrange that P(|Ci(G„)| = n) is 
only polynomially small in n. 

The next lemma is the combinatorial heart of the proof of Theorem 14.91 
Unfortunately, we cannot bound the expectation of x directly, only the con- 
tribution from components up to size some small constant times n. Formally, 
given a graph G with n vertices and a (5 > 0, let 

UG):=- nv)\ = - E IC^I'- (4-12) 

v<iy {G) ■.\C{v)\<in i:\Ci\<&n 

Note that if IC2I < (5n < |Ci|, then Xi^G) = x{G). 

Given a kernel k and an M > 0, we write k^^ for the pointwise minimum 
of K and M. 

Lemma 4.13. Let k be an irreducible, integrable kernel on an i.i.d. vertex 
space V with ||Tk|| > 1, and let s > and M > 0. Then there is a 6 = 
6{e, M, k) > such that 

00 

Ex5(G^(n, k)) < /i(5)-i ^(7^1, 1)^ + o(l), 

j=0 

where fi is the measure on S defined by djX{x) = f{x) dfi with 

f[x) = (1 - p((l - e)K*^ x) + 5e) A 1, (4.13) 
and is the integral operator on {S,fi) with kernel k. 
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Proof. As usual, we may and shall assume that = 1. 

Note that the statement becomes stronger if we increase M and/or de- 
crease e. Thus we may assume that (1 — e)K^'^ is supercritical, and that 
p((l - e)K^O > 2e. We also assume that M > 1 and e^^ < 1 + 5e. 

Let < 5 < e/M be a small constant to be chosen later, depending only 
on K, e and M, and let = nxsiGn) denote the number of ordered pairs 
(vjw) of vertices of Gn = G^{n,K) such that v and w are in a common 
component of size at most 5n. Also, let Nj denote the number of such pairs 
joined by a path of length j. Since < Yl^j^^ suffices to show that 

for < j < (5n we have 

EiV,-/n< (l|'l,l)^ + o(l/n), (4.14) 

with the error bound uniform in j. 

We may bound Nj by the number of paths of length j in G„, lying in 
components with at most 5n vertices. Thus '^Nj is at most n^~^^ times the 
probability that 12 • • • (j + 1) forms such a path. Let V' consist of the last 
(1 - e/M)n vertices of On- Coupling Gn and = G^(n, k^^) in the usual 
way so that G^ C let G' be the subgraph of G^ induced by V , noting 
that G' C Gn- Let A = Aj be the event that 12 • • • (j + 1) forms a path in 
Gn, and let B = Bj be the event that some vertex in [j + 1] is joined by an 
edge of G*^ to some component of G' of order at least 6n. Then 

ENj < n^+^F{AnB^). 

Unfortunately, we cannot quite prove the estimate we need for the right 
hand side above, so instead we use the less natural but stronger bound 

ENj< (^."^^EiN'^lBc), (4.15) 

where Nj is the number of ordered pairs {v, w) of vertices in Vq = [j + 1] 
such that V and w are joined in Gn by a path of length j lying in Vq (and 
thus visiting all vertices of Vq). 

Roughly speaking, the idea is to show that with very high probability 
Ci(G') will contain almost the 'right' number of vertices of each type, so that 
given the type y of one of the first j + 1 vertices, its probability of sending 
an edge to Ci(G') is almost what it should be, namely p((l — e/M)K*^;y). 
Unfortunately we cannot achieve this for all y, but we can achieve it for 
{xi, . . . , Xj+i}, which is all we need. Also, rather than working with Ci(G'), 
we work with the union of all components of order at least 6n. 

Let n' = (1 — e/M)n. Ignoring the irrelevant rounding to integers, G' has 
the distribution of G^{n\ (l-e/M)^*^), which dominates that of G"'^(n', (1- 

Recall that {\—£)k^^ is supercritical and that p{{l—e)n^'^) > 2e. Applying 
Lemma [412] to G' = G^(n',(l - £)k^^) we find that there is some c > 
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such that for any measurable A C S we have 

¥{Ci{A;G') < {n'{A)-2e/M)n) < F{Ci{A;G') < (n' {A)-e/M)n') < e'^^", 

(4.16) 

where ^' = ^(i_^)^a/. 
Let 

60 = min{e/M, 1/10} > 0, 

and fix < 5 < 5o chosen small enough that 

{e/SY < e^/\ (4.17) 

Let L denote the union of all components of G' of order at least 6n, and let 
L{A) be the number of vertices in L with types in A. If n'{A) > Ss/M and 
Gi{A; G') > {p'{A) — 2e/M)n, then since the final quantity is at least 5n we 
have L{A) > Ci{A; G'). Using (ITOD . it follows that 

F{L{A) < {fi'{A) - 3e/M)n) < e"^" (4.18) 

for any A; the condition is vacuous if f^'{A) < 3e/M. 

Given y £ S and i > 0, let Ay^i = {x G 5 : k^'^ {x, y) > ei}. Let £y be the 
event that L{Ay^i)/n > n'{Ay^i) - 3e/M holds for all i with 1 < i < M/e. 
Applying (j4.18p M/e = 0(1) times, we see that 

P(£:^) < (M/e)e-^" = 0(e"^"). (4.19) 

If £y holds, then 

M/e M/e M/e 

^K*^(x„,y) > ^L(^y,i)e > ^e(/i'(Aj,,i)-3e/M)n > n ^ e/i'(Aj,,i)-3en. 

v(^L i=l i=l 1=1 

Now Ay^i is empty for i > M/e, so we have 

M/e 00 „ 

Te^l'{Ay,,) = Te^i'{x■.K^'{x,y)>e^}= / e[K^' {x,y)/e\ dfx' (x) 

> / K^^{x,y)dfi'{x)-e= [ K'^{x,y)p{{l-e)K^;x)dfi{x) -e. 
Js Js 

Putting these bounds together, writing k' for (1 — e)K*^, we have 
^K^^(y,x„)/n > j K^'^ {x,y)p{K;x)dfi{x) - Ae 

= {T,Mp^,){y) - 4e > {T^>p^'){y) - 4e. 



Recalling that k! is supercritical, from (j2.14p we have Tf^^'P^' = — log(l — /)«;'); 



so when £y holds we have 



^^'{y. ^v)/n > - log(l - p{k'; y)) - 4e, 
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and hence 

l[il-K^\y,x,)/n) < {l-p{K';y))e^' < i _ y) + 5^. 

Since is bounded by M, and the product is always at most 1, it fohows 
that if £y holds and n > M, then 

ll{l-{K'^\y,x,)/nAl))<f{y). (4.20) 

Let £ = £xi n • • • n Sxj-f-i- Note that G' is independent of xi, . . . ,Xj+i. 
Given these types, from (j4.19p we have F{£) = l — 0{je~'^^) = l — 0{ne~'^"'), 
with the implicit constant independent of the types. Hence, we have F{£) = 
1 — 0{ne~'^'^) unconditionally. Then, for j < 6n, 



n 
.j + 



(4.21) 

using (j4.17p in the last step. 

Estimating N'- by the number of paths of length j lying in Vq, 

lE(A^jlBm^) < (^^. ^^^(j + l)!P(^nB'=nf) <n^+^¥{AnB^n£). 

(4.22) 

To estimate the final probability let us condition on G' and also on the vertex 
types xi, . . . assuming as we may that £ holds. Note that we have 

not yet 'looked at' edges within Vq, or edges from Vq to V . The conditional 
probability of A is then exactly 

i j 
Y\_{K{xi,Xi+i)/n A 1) < JJ K{xi,Xi+i). 

i=l 1=1 

For each i < j + 1, since f^. holds we have from (j4.20p that the probability 
that i sends no edge to L is at most f{xi). These events are (conditionally) 
independent for different i, so 

j i+i 
F{AnB''n£ \ xi,..., Xj+i) < J] ^(xi, x^+i) J] f{xi). 

i=l i=l 

Integrating out we find that 

. j i+i 
n^+i ¥{Af^B^^^£) < n W K{xuxi+i) JJ f{xi) d^i{xi) ■ ■ ■ d^i{x,+i) 

Jsi+^ i=l i=i 

= ri{Tll,l)^. 

From (021) it follows that (^"J E(iVjlf5cn£:) < 1)^. Combined with 

(j4.2ip and M.lSh this establishes (j4.14p : as noted earlier, the result follows. 

□ 
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Taking, say, M = 1/e and defining /e(x) by (14.130 . as e ^ we have 
(1 — e)K*^ y K pointwise, and hence p{{l — e)«;*'^;x) p{k;x) pointwise. 
Thus fe{x) \ 1 — p{k; x) pointwise. If we know that (T^l, < oo for some 
e > 0, then by dominated convergence it follows that (T^l, l)/i \ (T^l; 
Furthermore, if we have 



oo 



^{Til,l)f,<oo (4.23) 

j=0 

for some e > 0, then by dominated convergence, as e — > we have 

oo oo 
j=0 j=0 

Unfortunately we need some assumption on k. to establish ()4.23p . 

Proof of Theorem \4-9[ Suppose for the moment that (j4.23p holds for some 
e > 0, where jl is defined using fe{x), which is in turn given by (j4.13p with 
M = say. 

By the comments above, it follows that, given any ?? > 0, choosing e small 
enough and M large enough we have Yl'jLoi'^i^^ — x(^)+^- Lemma [4.13l 
then gives E xsiGn) < x(^) + 2rj if n is large enough, for some 5 = 5{r]) > 0. 
Hence, if 6 = 5{n) tends to zero, we have 

limsupEx5(G'„) < H'^). (4.24) 

Since k is supercritical we have p(k) > 0, and by (I2.2ip we have |Ci(G„)| > 
p{K)n/2 whp. For any fixed (5 > 0, by (12.220 we have |C2(Gn)| < 5n whp; this 
also holds if6 = 5{n) tends to zero sufficiently slowly. Given a function 6{n), 
let £n be the event that |C2(Gn)| < n5{n) < |Ci(G„)|. Then, provided 5{n) 
tends to zero slowly enough, £"„ holds whp. When £"„ holds we have XsiGn) = 
X(G„), soEixiGn);£n) < ExsiGn), and ^M) gives limsupE(x(G„); £:„) < 
x('t)- By Lemma 14.41 this implies that xiGn) — ^ x('^)) which is our goal. 
It thus suffices to establish that (|4.23p holds for some e > 0. 

Recall that /e(x) < 1 and \ /o = 1 — as e — > 0. Recall also that 
is defined as the integral operator 

9^ j i<-ix,y)giy)dp{y) = j K,{x,y)fs{y)g{y)dfi{y) 

on L?'{p). The map g{x) ^ g{x)fs{x)^^'^ is an isometry of L?{p) onto L^(/i), 
and thus is unitarily equivalent to the integral operator on Lp'{p) 
with kernel f^ix)^^"^ K,{x,y)fs{yy^'^ ■ In particular, ||Tk|| = ||re||, and for the 
special case e = 0, when = Tj^, llr^H = ||To||. 

Fix 6 > 0. Since is compact, there is a finite rank operator F with 
II A II < 5, where A = — -F. Let and A^ denote the operators obtained 
by multiplying the kernels of F and A by fs{x)^^'^fs{y)^^'^- Since /e < 1 
holds pointwise, we have 

IIAell < ||A|| < 5. 
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For any g £ the pointwise product f^g converges to fog in L^. Since F 
has finite rank, it follows that — FqII — *■ 0, and hence that 

limsup llTe — Toll < limsup — Fq\\ + 6 = 6. 

Since 5 > was arbitrary, we have ||T£ — To|| — > 0, and in particular ||Tk|| = 
||Te|| — > ||To|| = ||Tk|| < 1. Hence, there exists e > such that ||Tk|| < 1. 
But then holds, because (t|i, 1)^ < ||T«,|p'. □ 

Remark 4.14. Chayes and Smith [3] have recently proved a result related 
to Theorem 14. 7( i) or Theorem 14.81 foi' the special case where the type space 
S is finite. Their model has a fixed number of vertices of each type, which 
makes essentially no difference in this finite-type case. Chayes and Smith 
consider (in effect) the number of ordered pairs {vjw) of vertices with v 
of type i, w of type j, and v and w in the same component, normalized 
by dividing by n, showing convergence to the relevant branching process 
quantity. These numbers sum to give the susceptibility, so such a result is 
more refined than the corresponding result for the susceptibility itself. 

In our setting, the analogue is to fix arbitrary measurable subsets S and 
T of the type space, and consider xs,TiGn), which is 1/n times the number 
of pairs (f , w) in the same component with the type of v lying in 5* and that 
of w in T. The corresponding branching process quantity is just xs,t{h), 
i.e., the integral over x £ S of the expected number of particles in Xfi{x) 
with types in T. In analogy with Theorem 13.31 in the subcritical case this 
quantity may be expressed as X5,t('^) = {{I — Tk)~'^1s, It)/^ < oo. It is not 
hard to see that the proof of Theorem 14.81 in fact shows that 

XS,T{Gn) ^ X5,t(«), (4.25) 

where Gn = G^(n, k) is defined on an i.i.d. vertex space. The key point 
is that, in the light of Theorem 14.11 and its proof, it suffices to prove a 
convergence result for the contribution to XS,t(G„) from components of a 
fixed size k. For all the models we consider here, this may be proved by 
adapting the methods used to prove convergence of Nj.[Gn)/n; we omit the 
details. Once we have such convergence, we also obtain the analogue of 
()4.25p for so all our results in this section may be extended in this way, 
with theproviso that when considering G^(n, k) with a general vertex space 
V as in [J|, we must assume that S and T are ^-continuity sets. 

Remark 4.15. We believe that all the results in this section extend, with 
suitable modifications, to the random graphs with clustering introduced 
by Bollobas, Janson and Riordan [5j, and generalized (to a form analogous 
to G{An)) in 01; these may be seen as the simple graphs obtained from an 
appropriate random hypergraph by replacing each hyperedge by a complete 
graph on its vertex set. Note that in this case the appropriate limiting object 
is a hyper kernel (for the defintions see [^), and the corresponding branching 
process is now a (multi-type, of course) compound Poisson one. 
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A key observation is that in such a graph, which is the union of cer- 
tain complete graphs, two vertices are in the same component if and only if 
they are joined by a path which uses at most one edge from each of these 
complete graphs. Roughly speaking, this means that we need consider only 
the individual edge probabilities, and not their correlations, and then ar- 
guments such as the proof of Theorem 14.81 and (at least the first part of) 
Theorem 14. 71 go through with little change. It also tells us that the suscepti- 
bility of a hyperkernel is simply that of the corresponding edge kernel; this 
is no surprise, since for the expected total size of the branching process all 
that matters is (informally) the expected number of type y children of each 
type X individual, not the details of the distribution. This does not extend 
to the modified susceptibility Xi since this depends on the (type-dependent) 
survival probability /3(k; x), which certainly is sensitive to the details of the 
offspring distribution. 

Adapting the proof of Theorem 14.91 needs more work, but we believe it 
should be possible. Most of the time, one can work with bounded hyper- 
kernels, where not only are the individual (hyper)matrix entries uniformly 
bounded, but there is a maximum edge cardinality. Taking the r-uniform 
case for simplicity, one needs to show that the number of (r — l)-tuples of 
vertices in the giant component in some subset of is typically close to 
what it should be, since, in the proof of Lemma 14.131 t^is sets Ay^i should 
(presumably) be replaced by corresponding subsets of S^~^. For strong 
concentration, one argues as here but using the appropriate stability result 
from 0] in place of Theorem l4.11[ Needless to say, since we have not checked 
the details, there is always the possibility of unseen complications! 

5. Behaviour near the threshold 

In this section we consider the behaviour of x s-^icl x fo'^ ^ family \k of 
kernels, with n fixed and A ranging from to co. Since IITxkII = A||Tk||, then, 
as discussed in |4l], \k is subcritical, critical and supercritical for A < Acr, 
A = Acr and A > Acr, respectively, where Acr = ||^ft||~^- Note that if \\Ti^\\ < 
oo, then Acr > 0, while if HTkH = oo, then Acr = 0, so Xk is supercritical for 
any A > 0. 

Note also that Theorem 13 . 5 1 provides an alternative way of finding Acr (and 
thus IITkII): we can try to solve the integral equation / = 1+Tak/ = l+AT^/ 
and see whether there exists any integrable positive solution. This tells us 
whether xi^^^) is finite; since (by Theorems 13.31 and 13. 4p the susceptibility 
is finite in the subcritical case and infinite in the supercritical case, this 
information determines Acr. The advantage of this approach over attempting 
to solve (I2.14P itself is that the equation is linear; this is one of the main 
motivations for studying x- (Another is that it tends to evolve very simply 
in time in suitably parameterized models.) 

In the subcritical case, A < Acr, we have the following simple result. 
(When we say that a function / defined on the reals is analytic at a point x, 
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we mean that there is a neighbourhood of x in which / is given by the sum 
of a convergent power series; equivalently, / extends to a complex analytic 
function in a complex neighbourhood of x.) 

Theorem 5.1. Let k be a kernel. Then A ^ x{^i^) = xC-^^) 0,^ increas- 
ing, analytic function on (0, Acr), with a singularity at Acr- Furthermore, 
x(Ak) / x(AcrK) = x(AcrK) < oo as \ / Acr, and x(Ak;j;) / xiKr^^x) 
pointwise. 

Proof. By 

oo 

x(Ak)=M5)-^ j;(Til,l)A^ (5.1) 

j=0 

which converges for < A < A^ by Theorem l3.3[ Hence, xi^^^) is increasing 
and analytic on (0, Acr)- Moreover, by Theorem I3.4( ii). the sum in ()5.ip 
diverges for A > Acr; hence the radius of convergence of this power series is 
Acr- Since the coefficients are non-negative, this implies that x('^'^^) is iiot 
analytic at Acr- 

Finally, x(-^^) /" x('^cr'^) as A ^ Acr by (15. Ih and monotone convergence. 
Similarly, xi^i'^] ^ xl-^cr^! x) by (13. Ih and monotone convergence. □ 

We shall see in Section [Ql that it is possible to have x(Acr/^) < oo. As 
we shall now show, if is compact, then x('^cr^) = oo, and the critical 
exponent of x is — 1, as A Acr. 

Theorem 5.2. Suppose that is compact (for example, that f < oo). 
Then for some constant a, < a < 1, we have 

x{Xk) = x(A/^) = + 0(1), < A < Acr, 

and x(AcrK) = xi^cru) = oo. 

If, in addition, k is irreducible, then a = {J^ip) / where tp is any 

non-negative eigenfunction of T^. 

Proof. Since a compact operator is bounded, Acr > 0. We may assume that 
li{S) = 1 by Remark 12.11 Furthermore, we may replace k by Acr^ and may 
thus assume, for convenience, that ||Tk|| = 1 and Acr = 1- 

Let El be the eigenspace {/ G ^^(a*) : ^k/ = /} of T^, and Pi the 
orthogonal projection onto Ei. Since is compact and self-adjoint, Ei and 
its orthogonal complement are invariant, 1 does not belong to the spectrum 
of T^ restricted to E^ , and, for A < 1, ||(/ - \T^)-^{I - Pi)\\ = 0(1), while 
(/ - XT^y^Pi = (1 - A)-iPi. Consequently, by Theorem [231 

x(Ak) = (1-A)-1(Pi1,1)+0(1). 

Let a := (Pil,l) = HAlHi > 0; then a < = 1, so < a < 1. If 

a = 0, then Pil = 0, so the constant function 1 is orthogonal to Ei. But 
this contradicts the fact that Ei always contains a non-zero eigenfunction 
> 0, see the proof of Theorem 13.41 and [1, Lemma 5.15]. Hence, a > 0. 
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The fact that xi^crf^) = c« now follows from Theorem 15.11 
Furthermore, if k is irreducible, then Ei is one-dimensional, see again [1, 
Lemma 5.15 and its proof], so Pif = ||V'||2 ^(/) V')''/'> and the formula for a 
follows, noting that every non-negative eigenfunction is a multiple of this 
jp. □ 

In the supercritical case, only x is of interest. If we allow reducible k, we 
can have several singularities, coming from different parts of the type space, 
see Example 16.81 We therefore assume that k is irreducible. Even in that 
case, it is possible that the dual kernel k is critical, see [3, Example 12.4]; 
in this example it is not hard to check that x(^) is infinite. 

We conjecture that when k is irreducible, xi^i^) is analytic for all A 7^ 
Acr under very weak conditions, but we have only been able to show this 
under the rather stringent condition (15.21) below. (See also the examples in 
Section [H) Under this condition, we can also show that the behaviour of 
X is symmetric at Acr to the first order: the asymptotic behaviour is the 
same at the subcritical and supercritical sides. As seen in Examples 16.21 
and 16.31 this does not hold for all k, even if we assume the Hilbert-Schmidt 
condition J < 00. (Furthermore, we shall see in Sections 16. II and [6?2] that 
the second order terms generally differ between the two sides.) 

Theorem 5.3. Suppose that k is irreducible, and that 

X Js 

(i) The function A ^ xi^f^) analytic except at A^ := HTkH"-^. 
(h) As \^ Acr, 

x(Ak) = -^ + 0(1), 

I A — Acr I 

with b = (/^ ^)^/ ip'^ > 0, where ip is any non-negative eigenfunc- 
tion ofT^. 

Proof. The subcritical case A < Acr follows from Theorem 1 5. 2 1 so we assume 
A > Acr- (Note that (j5.2p implies that is Hilbert-Schmidt and thus 
compact.) We may further assume that /i(5) = 1. 

(i): Let Aq > Acr- By |4, Section 15], there exists an analytic function 
z ^ p'^ defined in a complex neighbourhood U of Aq and with values in the 
Banach space L'^{p) such that = pzn when z is real, and (j2.14p extends 
to 

p+ = l-e~^^^PK (5-3) 

We may further (by shrinking U) assume that ||/5^||2 is bounded in U . Then, 
by (|5.2p and Cauchy-Schwartz, ||T'k(p+)||oo = 0(1) in [/, and thus, by (j5-3p . 
|1 — p^l is bounded above and below, uniformly for z £ U. In particular, for 
every Xk with real X £ U, L'^ip) = L'^{p), with uniformly equivalent norms- 
We can therefore regard as an operator in Lp'{p,). 
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We define, for zGU, fj := zT^{{l- p+)f); tlius fx = T^^ for real A G C/ 

by ()2.16p . Note that z is an analytic map of U into the Banach space 

of bounded operators on 

By Theorem 13.31 ^ ~ is invertible. By continuity, we may assume 

that I — Tz is invertible in U. Then f{z) := {{I — T^)"^!, 1 — p^)^ is an 
analytic function in U, and /(A) = xC-^/^) ^oi real X G U hy Theorem 13. 3f ii). 
Hence xi^f^) is analytic at Aq. 

(ii): We use a result from perturbation theory, for convenience stated as 
Lemma [531 below in a form adapted to our purposes; see [H, Section VII. 6] 
or [l^] for similar arguments and many related results. 

We may rescale and assume that Acr = = 1, i-e., n is critical. 

It will be convenient to use the fixed Hilbert space L'^ifJ-) rather than 
-L^(/2); recall that ju depends on A. Define a self-adjoint operator Tx in 
L\fi) by 

Txf := (1 - PA.)'/'Ar,(/(l - px.)'/'), (5.4) 
and note that if Ux is the unitary mapping / (1 — pxkY^"^ f of L^{p) 
onto -L^(/u), then Tx = UxT^U^^ by (I2.16p . Hence, Tx in -^^(a*) is unitarily 
equivalent to T-q in Lp'{p). Further, by Theorem 13. 3( ii). 

xiXt^) = {{I - T^r\ l)n = iil - Txr'Uxl, Uxl)^. (5.5) 

Note that p^ = 0, and thus = T^, which has a simple eigenvalue 1, 
with a positive eigenfunction ij; [J, Lemma 5.15], and all other eigenvalues 
strictly smaller. We may assume that HV'lb = 1- 

We apply Lemma 15.41 with T = Ti and T' = Tx, with A = 1 + e for small 
e > 0. By \i, Section 15], || /OAkIIoo — 0{e), and more precisely, pxk — OeV'+Pe 
with ||p*||2 = O(e^) and 

ae= r ,^3, g + 0(£'). (5.6) 
Js ^ 

It follows (recalling that ij: is bounded because ip = T,^ip and (j5.2p ) that 
(1 - pxnY^'^i' = ip - ^aeilP' + Tg, with ||r£||2 = O(e^). Consequently, (|5.4p 
implies that ||TA-fi|| = 0(e) and, using = = = 

V^d/i and (l5^ . 

(T^V, V') = A(r4(i - PA«)'/V) , (1 - Px.f'^i') 

= A((r,V,V') - \ae{T^^,i>^) - ia,(r«V'',V') +0(e2)) 
= {l + e){l-2e + 0{e^)) 
= l-e + 0{e^). 

Further, Uxl = (1 - P\kY'^ = 1 + 0(e)- Hence, ([53]) and ([521) yield 

(M:)2 + o(e) M)!^nn^ 
X((l + e)K) = ^^oie'^) ^ ' " — e — 
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which is the desired result. □ 

Lemma 5.4. Let T be a compact self-adjoint operator in a Hilhert space H , 
such that T has a largest eigenvalue 1 that is simple, with a corresponding 
normalized eigenvector tp. Then there exists rj > such that if T' is any 
self-adjoint operator with \\T' — T\\ < rj such that I — T' is invertible, then 

for any f,g £ H with \\f\\, \\g\\ < 1. 

Proof. The spectrum cr(T) C (— oo, 1 — 5] U {1} for some 5 > 0. Let 7 be the 
circle {z : \z — 1\ = 5/2}. Then, as is well known, the spectral projection 

Po := f^izl - T)-Uz (5.8) 

is the orthogonal projection onto the one-dimensional eigenspace spanned 
by ■0. Let A = T' — T. If ^ is any self-adjoint operator with \\A\\ < rj, for 
some sufficiently small r/ > 0, then zl — T — A is invertible for 2: G 7, and 
we define 

Thus Pa is the spectral projection for T -\- A associated to the interior of 7. 
It follows from ([52]) and dSj]) that \\Pa - PqW = 0{\\A\\), so if r] is smah 
enough, \\Pa — -Poll < and it follows fli, Lemma VII.6.7] that Pa too 
has rank 1; this must be the orthogonal projection onto a one-dimensional 
space spanned by an eigenfunction ipA oi T + A with eigenvalue Xa, with 
\Xa — 1| < d/2. Moreover, if Xa 7^ 1, then since all other eigenvalues of 
T + A then lie outside 7, 

(/ - (T + A))-^ = (1 - XaT^Pa + Ra, (5.10) 

with ||i?^|| < 2/6 = 0(1). 

Since ||PaV'-V'II = II(^'a--Po)V'II = O(PII), Pa^^ + (provided r? is smah 
enough), and thus we can take ■i/'A = -fUV'- Hence HV'^i — V'll = ll-fUV' ~ V'll = 
O(PII) and 

(Va,0) = (0,V') + O(P||) = 1 + O(P||), 

(rvA,0) = k^A.m = (V'A,V'> = i + o(Pll), 

{Ai;A,i^) = {Ai^,i^) + 0{\\Af), 

and thus 

A, = (JI±M±^1. = 1 + = 1 + (^^,,,,) + o(Un (5.11) 

The result follows from ()5.10p and ()5.1ip . using Pq/ = {f,'^)'^- D 
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6. Examples 

In this section we give several examples illustrating the results above and 
their limits. We sometimes drop k from the notation; we let pk denote the 
function pk{x) = pk{K;x). (But we continue to denote the number J^p^dp 
by Pki^), in order to distinguish it from the function pk-) 

Note first that the probabilities pk{x) can in principle be calculated by 
recursion and integration. The number of children of an individual of type 
x in the branching process is Poisson with mean J K{x,y) dp{y) = T^l{x), 
and thus (in somewhat informal language) 

pi{x) = P(x has no child) = e~^''^'^''\ (6.1) 

Next, |Xk(x)| = 2 if and only if x has a single child, which is childless. 
Hence, by conditioning on the offspring of x, 

p,(x) = e-^«i(-) / K{x,y)¥{\X,{y)\ = l)dp{y) = e-^«iWrjpi)(x) 

= pi{x)TM{x). (6.2) 
Similarly, considering the two ways to get |Xk(2;)| = 3, 

Pi{x) = e-^''^^''^ / K{x,y)p2{y)dp{y) 



+ e ^''^("^^ ]^ j^K{x, y)pi (y) dp{y) k{x, z)pi (z) dp{z) 

= pi{x)TM{x) + \pi{x){TM{x))\ (6.3) 
and the three ways to get |Xk(x)| = 4, 

PA = PiTips) + piT{pi)T{p2) + \pi{Tpif, (6.4) 

and so on. In general, for p^, k >2, we get one term pi Y\j T(pj)"^^ /rrijl for 
each partition l™i2'^2 . . . of fc - 1. 

The numbers Pkii^) are then obtained by integration. 

Alternatively, a similar recursion can be given for the probability that 
Xk(x) has the shape of a given tree; this can then be summed over all trees 
of a given size. 

6.1. The Erdos— Renyi case. Let S consist of a single point, with p{S) = 
1. Thus, K is a positive number. (More generally, a constant k on any 
probability space (5, p) yields the same results.) We keep to more traditional 
notation by letting k = A > 0; then G(n, k) = G{n,p) with p = X/n. See [1, 
Example 4.1]. 

Since is just multiplication by A, ||Tk|| = A, and, as is well-known, k 
is subcritical if A < 1, critical if A = 1, and supercritical if A > 1. 
In the subcritical case, by (13. 2p or Theorem 13. 3( i). 

x{k) = A < 1. (6.5) 
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Theorem 14.71 or Theorem 14.81 shows that x{G{n,\/n)) — ^ (1 — A)~^ for 
every constant A < 1. (This and more detailed results are shown by Janson 
and Luczak [g^] by another method. See also Durrett (itI . Section 2.2] for 
the expectation Ex(G(n, A/n)).) 

Similarly, if A > 1 then x(G(n, A/n)) — ^ x{^^) = cxd by Theorem 13.41 and 
any of Theorems 14.51 14.71 or 14.81 

For X, we have the same results for A < 1. In the supercritical case A > 1, 
Tk is multiplication by A(l — p(A)) < 1, where 1 — p{\) = exp(— A/9(A)) by 
(imil . Hence, by Theorems SiT] and EH or ([331), for A > 1, 

x(G(n, X/n)) x(-) = ^ _ ^(i _ ,(a)) = l-A(l-p(A)) - ^^'^^ 

More generally, Theorem 14. 71 shows that x(G(n, A„/n)) — ^ x(A) for every 
sequence A„ — > A > 0. 

For A = 1 + e, e > 0, we have the Taylor expansion 

8 2 28 , 464 4 , 

P(l + e)=2e-je^ + -e»--e* + ... (6.7) 

and thus 

X 1 4 4 176 , , , 

Combining (j6.5|) and (|6.8|) . we see that, as shown by Theorem 15.31 x(A) ~ 
1/1 A — 1| for A on both sides of 1, but the second order terms are different 
for A / 1 and A \ 1. 

We can also obtain x(A) and x(A) from pk and the formulae p.7p and 
()2.8p . In this case, is an ordinary, single- type, Galton- Watson process 
with Poisson distributed offspring, and it is well-known, see e.g. 27; 33; 



18|; [3^; 123] , that |Xk| has a Borel distribution (degenerate if A > 1), i.e., 

; k—1 

p^,{K) = pu{x) = -^\^-^e-'^\ k>l. (6.9) 
Consequently, if T(z) := YlV=i ^HtT'^^ ^^^^ function, then 

p{^) = l- Pk{n) = l-^^^^ (6.10) 

l<fc<oo 



and, using the well-known identity zT'{z) = T{z)/{\ — T{z)), see e.g. [2l| . 

l<fc<oo k=l ^ ' 

In the subcritical case, when A < 1, we have T{\e~'^) = A, and we recover 
(I63D . In general, (l6T0]) and (l6TT]) yield (fO]) . 

Remark 6.1. Consider the random graph G{n,m) with a given number 
m of edges. In the subcritical case m ~ An/2 with < A < 1, we obtain 
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x{G{n, m)) — ^ x{i^) = 1/(1 ~ W comparison with G{n,p) with p = Xn/n 
for A„ = 2m /n it n~^^^, say, using Lemma 12.21 In the supercritical case 
A > 1, one can use standard results on the numbers of vertices and edges in 
the giant component; conditioning on the giant component assuming typical 
values, the rest of the graph is essentially a subcritical instance of G{n, m) 
with different parameters; this may be compared with G{n,p) as above. 
Consequently, for m ~ An/2 with A > 1, x{G{n,m)) — ^ xi'^)^ where x(^) 
is given by (j6.6p and (|6.1ip . just as for G{n,p) with p = X/n. 

6.2. The rank 1 case. Suppose that k{x, y) = ip{x)ilj{y) for some positive 
integrable function ^|J on S. This is the rank 1 case studied in 0, Section 
16.4]; note that is the rank 1 operator / i-^ {f,ip)tl', with ip as eigenfunc- 
tion, provided ^p G -^^(a*)- 

We assume, for simplicity, that /u(5) = 1. As in Section[5]we consider the 
family of kernels Xk, A > 0. In this case, ||Tk|| = HV'lli = /^V'^) aiid thus 

Acr = IIV'lls 

In the subcritical case, A < A^ = (/ V'^) > which entails ifP' < oo, we 
have by induction 

Tll{x)=X^U V'd/i)'"' / Vd/.-V'(x), J>1, 



and thus by (j3.2p (or by solving (j3.5 



2 



x{Xk) = x{Xk) = 1 + , /,„ = 1 + 



l-A/V'2 1-A/Acr 

In particular, this verifies the formula in Theorem 15. 2[ 

In the supercritical case, we first note that the equation (j2.14p for p = pxn 
becomes 

p = I - e-^'^'^P = 1 - e-^^P^^^"^ . (6.13) 
We define ^ G (0, cxd) by ^ := X{p,ip), and thus have 

p = l- e-^'^, (6.14) 

with ^ given by the implicit equation 

C = X f p{x)ij{x) dp{x) = X j V'(a;)(^l -e-«^(^)) d/i(x). (6.15) 



(See [J, Section 16.4], where the notation is somewhat different.) We know, 
by results from [3], that (j6.13p has a unique positive solution p for every 
A > Acr; thus (j6.15p has a unique solution = ^(A) > for every A > A^. 
It is easier to use ^ as a parameter; by (j6.15p we have 

^ = /(l-e-e^)V'- ^'-''^ 
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The denominator is finite for every ^ > since tp & L^; moreover, /(I — 
e~^'^)'ip < J S.ip'^, and thus (|6.16p yields A > l/JV'^ = Acr- Consequently, 
(j6.15p and (j6.16p give a bijection between A G (Acr,c>o) and ^ G [0, oo). 
Furthermore, differentiation of ()6.16p shows that A = A(^) is differ entiable, 
and it follows easily from /(I - e~^^)ip > J ^-i/'^e"^'^ that dA/d^ > 0. 
Hence, the function A(.^) and its inverse ^(A) are both strictly increasing and 
continuous. In particular, A \ A^ <^=^ \ 0. Moreover, the denominator 
in (j6.16p is an analytic function of complex ^ with Re^ > 0; hence A(^) and 
its inverse ^(A) are analytic, for ^ > and A > Acr, respectively. 

We note also the following equivalent formula, provided J^ip"^ < oo: 

T^-T = r'/ fe-«^ - 1 + CV') V'- (6.17) 



Acr A 

By (I2J6D and (imi) . 

J s 

(6.18) 

Hence T-^ too is a rank 1 operator, with eigenfunction ^ and eigenvalue 
(take / = ip in (f08]) ) 

7 = A e-^^^^H{xf Mx) = (6.19) 

Since y'^e~'^ < y{l - e"^) for y > 0, it follows that < 7 < 1. (When 
J ip'^ < 00, this follows also from the general result P, Theorem 6.7], cf. 
Theorem 13.31 ) Hence I — is invertible (in, for example, -L^(/2)), and by 
Theorem 13. 3( ii). 

xiXK; x) = {l- p{x)){I - T^J-'l{x) = e-«^(-)(/ - T^J^'l{x). (6.20) 

Let us write g := (I - T^J'^l. Then, by 1 = " T^Jd = 9 - d^, 

with C = A e~^^ipg. Hence, g = 1 + aiid, using (j6.19p . 



C = X e-^^iPg = X / e-^'^V + AC / e'^^V = A / e'^'^V + Ct 
Hence, using (j6.16p and (j6.19p . 



C 



1-7 /(l-e-?'/')V'-^/e-«^V^' 
Finally, by (lOOD . 

X = [ x{X>^;x)dp{x)= [ e-^^g= [ e'^'^ + C / e"?'^^ 



+ J(l_e-5^(l+CV'))V.- ^ 
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We observe that (I6.21|) shows that x is an analytic function of G (0, oo), 
and thus of A G (Acr, oo). (So in the rank 1 case, at least, the condition (j5.2p 
is not required for Theorem 15. 3( i).) 

Next, suppose that J^ip^ < oo. In this case, we can differentiate twice 
under the integral signs in ()6.16p and ()6.2ip using dominated convergence 
(comparing with f^'ip'^), and taking Taylor expansions we see that as — > 
we have 

^ = ,r 12 iJr /3^ (,2^ = + hrr^ + ^(^) (6.22) 

U V' - 2^ J V''^ + ) 2 (J ^2) 



and 



where we used ()6.22p in the last step. 

Note that (|6.12p and (|6.23p show that the behaviour of x at the critical 
point Acr is symmetrical to the first order: 

- 'a/1-ii • 

at least when j ijj'^ < oo. (This is the same first order asymptotics as given 
by Theorem 15. SI fii). but note that the latter applies only when ■i/' is bounded, 
since ()5.2p fails otherwise.) The second order terms are different on the two 
sides of Acr, though: if J '0'^ < oo, then carrying the Taylor expansions above 
one step further leads to 

^~ A/Acr-l 1^3 + 3(^^3)2 

+ o(l), A\Acr, (6.25) 

in contrast to (j6.12p for A < Acr. 

To see what may happen if J^Tp^ = oo, we look at a few specific examples. 

Example 6.2. Let 2 < g < 3 and take 5 = [1, oo) with d^(x) = qx^'^^^ dx, 
and take 'ip{x) = x; note that < oo if and only if p < g; in particular 

J g ilP' < oo but J^ip"^ = oo. By (j6.17p . and standard integration by parts of 
Gamma integrals, as ^ ^ we have 

-| 1 fOO POO 

— --=rW {e-^^ -l+Cx)qx-'^dx = q^'^-^ {e'y - 1 + y)y~'^ dy 

Acr A Ji 

POO 

^^5-2 / (e-?;_i+y)y-9dy = g^5-2r(l-g), 
Jo 
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or A — Acr ~ '7r(l — q)\cr^'' '^- Similarly, by another integration by parts, 
y (l - e-^^{l + eV')) d/i = J (l - e'^'^il + ^x))qx-'^ dx 

/"OO /"CXD 

/ {l-e~yil + y))y-^dy^qC'-' (l _ e-?^(l + dy 
J£, Jo 

T{3-q)=q{q-2)e-'Til-q), 



q-1 
and thus by (fOT]) . 



9 ^(/^^)' (/^)'^c^ ... 

^ Q(g-2)e'^^ir(l-g) (g-2)(A-Acr)' ^ 

which still has power —1, but differs by a factor {q — 2)~^ from the subcritical 
asymptotics in (j6.12p and Theorem I5.2[ Hence, (j6.24p does not hold in 
general without assuming J^ip^ < oo. (Although this integral does not 
appear in the formula.) 

Example 6.3. We see in Example 16 . 21 that x is relatively large in the barely 
supercritical phase when is only a little more than square integrable. 
We can pursue this further by taking the same S and ip, and dfi{x) = 
c(logx + l)~'^x~^dx with q > 1 and a normalization constant c. Similar 
calculations using ()6.17p and ()6.23p (we omit the details) show that, letting 
c denote different positive constants (depending on g), as ^ — > we have 
A - Acr ~ c(log(l/0)-(''-i) and x ~ c(log(l/e))'?, and thus 

x(Ak) ~ c(A - Aer)-''/('^-^\ A\AcD 

with an exponent —q/{q — 1), which can be any real number in (— oo, —1). 

Taking instead d^(x) = c(loglogx)~^(logx)~^x~^ dx, x > 3, we similarly 
find A — Acr ~ c(log log(l/^))^^ and x ~ c(log(l/^))(log log(l/^))^, and thus 

x(Ak) = exp I - ^ _ ^ I, A\Acr, 

with an even more dramatic singularity. Of course, this sequence of examples 
can be continued to yield towers of exponents. 

6.3. The CHKNS model. Consider the family of kernels Ak, A > 0, with 

K{x,y):=-^-l (6.26) 
X V y 

on 5 = (0, 1] with Lebesgue measure /x. We thus have 

TxJ{x) = x(^^-l^ £ f{y) dy + A (i - l) /(y) dy 

= - r f{y) dy + X [' ^dy-X [' f{y) dy. (6.27) 

^ Jo Jx y Jo 
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Remark 6.4. Equivalently, by a change of variable, we could consider the 
kernel A(e^'^^ — 1) on 5 = [0, oo) with dfi = dx; we leave it to the reader 
to reformulate results in this setting. 

This kernel arises in connection with the CHKNS model of a random 
graph introduced by Callaway, Hopcroft, Kleinberg, Newman and Strogatz 
[12l |. This graph grows from a single vertex; vertices are added one by one, 
and after each vertex is added, an edge is added with probability 5 £ (0, 1); 
the endpoints are chosen uniformly among all existing vertices. Following 
Durrett [l6l : 17], we consider a modification where at each step a Poisson 



Po((5) number of edges are added to the graph, again with endpoints chosen 
uniformly at random. As discussed in detail in [4, Section 16.3], this yields 
a random graph of the type G^{n,Kn) for a graphical sequence of kernels 
with limit Xk, where A = 25, on a suitable vertex space V (with S and 
/i as above). 

Let us begin by solving (13. Sp . If / = Tak/ + 1, then (|6.27p implies first 
that / E C(0, 1) and then / G C^O, 1). Hence we can differentiate and find, 
using (I6.27[) again, that 

fix) = {TxJYix) = -4 r fiy) dy. (6.28) 

X Jo 

With F{x) := f{y) dy, this yields F"{x) = -AF(x)/x^ with the solution 
F{x) = Cix"+ + C2X"-, where a± are the roots of a{a — 1) = —A, i.e., 

a-t = ^ ± \J\ — X; if A = 1/4 we have a double root a+ = a_ = 1/2 

and the solution is F{x) = Cix^/^ + C2x'^/'^\ogx. Hence any integrable 
solution of (|3.5p must be of the form /(x) = Cj^x'^+^'^ + C-X°'-~^ , or /(x) = 
C+x~^/^ + C_x~^/^logx if A = 1/4. Any such / satisfies (|6.28|) . and since 
l!^?Ff\ yields Tx^,f{l) = 0, it solves ([SSD if and only if /(I) = 1, i.e., if 
C+ + C7_ = 1 (C+ = 1 if A = 1/4). 

If < A < 1/4, then < a_ < 1/2 < a+ < 1, so the solution /(x) = 
x"+~^ is in L^(0, 1) and non-negative; by Corollarv 13.61 this is the unique 
non- negative solution in L? , and 

X /"^ « 1 1 2 1 - VI -4A 

x(Xk)= x°+-Mx = — = , . (6.29) 

' Jo a+ i + Vr^riA 2A ^ ^ 

(If we are lucky, or with hindsight, we may observe directly that x"+~^ 
is a solution of (|3.5p by (j6.3ip below, and apply Corollary 13.61 directly, 
eliminating most of the analysis above.) 

For A < 1/4, we have shown that xi^f^) is finite, so Xk is subcritical; thus 
Acr > 1/4. Since the right-hand side in (I6.29P has a singularity at A = 1/4, 
Theorem 15.11 shows that Acr > 1/4 is impossible, so we conclude that Acr = 
1/4. (Equivalently, ]]Tfj]] = 4.) This critical value for the CHKNS model has 
earlier been found by Callaway, Hopcroft, Kleinberg, Newman and Strogatz 
[l3 ] by a non-rigorous method, also using ()6.29p which they found in a 
different way; another non-rigorous proof was given by Dorogovtsev, Mendes 
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and Samukhin ij] , and the first rigorousproof was given by Durrett la; ll7 1 . 
See also Bollobas, Janson and Riordan where different methods were 

used not involving the susceptibility. The argument above seems to be new. 

By Theorem 15. 1^ we can let A Acr in (|6.29p , and see that the equation 
holds for A = Acr = 1/4 too; i.e., xl-^cr'^) = 2. 

We see also that in the (sub)critical case A < 1/4, x(Ak;x) = x'^+~^ . 

We have no need for the other solutions of (13. 5p , but note that our analysis 
shows that for A < Acr, the other non- negative, integrable solutions of (j3.5p 
are given by + C(x°~~^ — with C > 0. Similarly, although 

we have no need for the solutions of (jS.Sp for A > Acr, let us note that 
for the critical case A = Acr, the argument above shows that there is a 
minimal non-negative solution which belongs to but not to l?\ 

there are further solutions x~^/^ — Cx""*^/^ log x, C > 0. For A > 1/4, 
the roots a± are complex, and the only real integrable solution to (j3.5p is 
i(a;"+-^+x"--^) = Rex"+-^ = ^-^/^ cos((A- i)^/^ log x) , which oscillates; 
thus there is no finite non-negative solution at all. 

Before proceeding to x i^i the supercritical case, let us calculate p\, for 
small k. We begin by observing, from (j6.27p . that T\t^\[x) = — Alogx. 
Hence (16. ip yields 

/)i(Ak;x) =6^1°*^^ = x\ (6.30) 
Further, by (j6.27p . for every non-zero 7 > — 1, 



r..(x-) = ;^^(l-x-). (6.31) 



Hence (|6.2p yields 



P2{\k- x) = x^rA,(x^) = -^(x^ - x^^). (6.32) 

i + A 

Similarly, ()6.3p and ()6.4p yield 

, (2 + 3A)x3^ - 4(1 + 2A)x2^ + (2 + 5A)x^ ,„ 

P%{\k\x) = —77- — , (6.33) 

' 2(1 A)2(l -h2A) ^ ' 

and a formula for P4^{Xk; x) that we omit, and so on. By integration we then 
obtain 

= (i+A)^a+2A) ' ^'-''^ 

3A^ 

^^^^"^ = (1 + A)3(1 + 2A)(1 + 3A) ' ^'-''^ 
2X^(7 + 15X) 

P'^^""^ - (1 + A)4(1 + 2A)2(1 + 3A)(1 + 4A)- ^^'^^^ 
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It is obvious that each pk{XK;x) is a polynomial in with coefficients 
that are rational functions in A, with only factors 1 + jX, j = 1, . . . , /c in the 
denominator. Hence, each p{Xk) is a rational function of the same type. 

There is no obvious general formula for the numbers pk^Xn), but, sur- 
prisingly, they satisfy a simple quadratic recursion, given in the following 
theorem. This recursion was found by Callaway, Hopcroft, Kleinberg, New- 
man and Strogatz [12], using their recursive construction of the graph, see 
also 0, Chapter 7.11. (The Argument in is non-rigorous, but as pointed 



out by Durrett it is not hard to make it rigorous.) We give here 

a proof that instead uses the branching process, which gives more detailed 
information about the distribution of the 'locations' of the components. 

Theorem 6.5. For the CHKNS kernel (j6.26p . pkiXn) satisfies the recursion 

kX ^'^ 

PkiXn) = ^ ^pk-jiXK)pj{XK), k>2, (6.38) 
i=i 

with pi{Xk) = 1/(1 + A). Hence, for each k > 1, pi.{Xk) is a rational function 
of X, with poles only at j = 1 . . . ,k. 

Moreover, each function pk{x) = pk{XK; x) is a polynomial in x'^, with co- 
efficients that are rational functions of X, which can he calculated recursively 
by 

d 

x—pk{XK;x) = kXpk{XK;x) - jXpk-j{XK)pj (Xk; x) , k > 1, (6.39) 

i=i 

together with the boundary conditions pi{Xk;1) = 1 and pk{XK;l) = 0, 
k>2. 

Proof. Fix A > 0. To simplify the notation, throughout this proof we write 
K for the kernel so far denoted Xk. Let e £ (0, 1/2), say, and let X'^ be 
with all points scaled by the factor (1 — e); this is the branching process 
defined by S' := (0,1 - e], dp' := (1 - e)-^ dx and K'{x,y) := X{^ - l). 
In X^, the offspring process of an individual of type x has intensity 

K'{x,y)dp'{y) = X(—^ T~) = f^i^ ^ v) - -r~ ^ 2/ < 

VxVyl — £/ 1 — e 

This is less than the intensity in X^- We let K'{x,y) = if x > 1 — e or 
y > 1 — e, and define K"{x,y) = K{x,y) — K'{x,y) > 0. More precisely, for 
< X < 1 - e and < y < 1, 

Thus Xk(x) and X'^(x) may be coupled in the natural way so that X'^{x) C 
Xfi{x) in the sense that an individual in X'i^{x), of type z say, also belongs to 
Xk{x), and its children in Xk{x) are its children in X'^{x) plus some children 
born according to an independent Poisson process with intensity k"{z, y) dy; 



Pk-j{K- y) dy+0{e) = pk-j{K)+0{e). (6.42) 
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we call the latter children (if any) adopted. An adopted child of type y gets 
children and further descendants according to a copy of X^iy), independent 
of everything else. Note that this adoption intensity K"{x,y) is independent 
of x S S', and that the total adoption intensity is Jq k" {x, y) dy = e\+0{e^). 

Fix A; > 1. If [^^(a;)! = k, then either |X'^(x)| = k and there are no adop- 
tions, or |X^(a;)| = j for some j < k and there are one or more adoptions, 
with a total family size of k — j. If |X^(x)| = k, then the probability of some 
adoption is keX + O(e^), and thus 

P(|X«(x)| = k I \X'^{x)\ = k) = I - kXe + O(e^). (6.41) 

Now, suppose that |j£'^(3;)| = j < k. The probability of two or more adop- 
tions is O(e^). Suppose that there is a single adoption. If the adopted child 
has type y, the probability that this leads to an adopted branch of size 
k — j, and thus to |Xk(x)| = k, is pk-j{K] y). By (j6.40p . the adoption inten- 
sity K"{z,y) is independent of z as remarked above, and is almost uniform 
on (0, 1]; it follows that the probability that |Xre(x)| = k, given |X'^(a;)| = j 
and that there is a single adoption, by some individual of type z in X'^(x), 
equals 

Jo i^"{z,y)Pk-j{K;y)dy _ 
!oi^"{z,y)dy Jo 
Since the probability of an adoption at all is jeX -|- O(e^), we obtain 

P(|X,(x)| = k I |X;(x)| = j) = jXpk-j{K)e + 0{e^). (6.43) 

Consequently, for every k > 1 and x £ (0, 1 — e], 

fc-1 

Pk{K;x) = {l-kXe)pkiK;x) + '^jXpk^j{K)pj{K;x)e + 0{e'^). (6.44) 

i=i 

(The implicit constant in O here and below may depend on k but not on x 

or e.) Replace x by {l — e)x and observe that, by definition, |X'^((1 — e)x)| = 
|Xk(x)| and thus Pj{K'; (1 — e)x) = pj{K;x). This yields 

fc-1 

Pk{K; {I - e)x) = {l-kXe)pk{K;x) + '^jXpk-j{K)pj{K;x)e + 0{e'^). (6.45) 

i=i 

Letting e \ we see first that pk{K; x) is Lipschitz continuous in (0, 1), and 
then that it is differentiable with 

d ''^^ 
x-^Pk{i^]x) = kXpk{K;x) -^jXpk-j{K)pj{K,]x), k>l, (6.46) 

i=i 

which is (|6.39p in the present notation. 

For /c = 1, (j6.46p gives pi{k;x) = Cx^, for some constant C. For x = 1 
we have K{l,y) = 0, so the branching process Xk(x) dies immediately, and 
Pi(k; x) = 1. Thus Pi{k; x) = x^ as shown in (I6.30p . For k > 2, we note that 
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X — pk{K; X) = [xpk{K] X)] Q - / Pk{K] X) dx = - Pk{K). 







xpk{K; x) ^ as X — > or X ^ 1, because pk{K; x) < 1 — pi{k; x) = 1 — x^, 
and thus, integrating by parts, 

dx' 

Hence, integration of (j6.46p yields the recursion formula 

k-l 

(1 + kX)pk{K) = Y,j^Pk-M)p3i^^)^ k>2. (6.47) 
i=i 

Replacing j hy k — j in the right-hand side of (j6.47p and summing the 
two equations, we find that 

fc-i 

2{l + kX)pk{K) = J2{j + k-j)Xpk-j{K)pj{K), k>2, (6.48) 
i=i 

which is dlJSl). □ 

The susceptibility x was calculated for all A by Callaway, Hopcroft, Klein- 
berg, Newman and Strogatz using the recursion formula (|6.38p . see also 
Durrett [l^: [l7l|. We repeat their argument for completeness. 

Let G{z) := ^kLi Pk{^i^)z^ be the probability generating function of 
IXakI, defined at least for \z\ < 1. Note that in the supercritical case, |Xak| is 
a defective random variable which may be oo; we have G(l) = 1 — P(|Xak| = 
oo) = 1 - p{Xk). Further, G'(l) = x(Ak) < oo. 

The recursion (j6.38p yields, most easily from the version (j6.47p . 

G{z) + XzG'{z) = XzG'{z)G{z) + {1 + X)pi{Xk)z = XzG'{z)G{z) + z, (6.49) 
and thus 

'''<^) = x^(i%y W < ' 

In the supercritical case, G{1) < 1, and we can let z y 1 in (|6.50p . yielding 
x(Ak) = G'{1) = 1/A. (In the subcritical case, rHopital's rule, or differ- 
entiation of (j6.49p . yields a quadratic equation for G'(l), with (j6.29p as a 
solution; this is the method by which (|6.29p was found in [l^.) 

Summarizing, we have rigorously verified the explicit formula by Callaway, 
Hopcroft, Kleinberg, Newman and Strogatz (l^ : 



xiXn) = 2A ' " - 4' (6.51) 




Note that there is a singularity at A = 1/4 with a finite jump from 2 to 4, 
with infinite derivative on the left side and finite derivative on the right side. 
It is striking that there is a simple explicit formula for xi^''^) = while 
no formula is known for G{1) = 1 — p{Xk,). This is presumably related to the 
fact that x(^'^) may be found by solving the linear equation ()3.5p . whereas 
p{Xk) is related to the non-linear equation (12.141) . As A = l/4-|-e \ 1/4, 
p{Xk) approaches extremely rapidly, as exp(— (7r/2\/2)£~^''^ -l-O(loge)) 
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0]; the behaviour at the singularity is thus very different for G(l) and 
G'{1). 

Note also that, by (12. Sh . the discontinuous function x(^'^) is the pointwise 
sum of the analytic functions kpk^K). 

Remark 6.6. We can obtain higher moments of the distribution {pk{^i^))k>i 
of by repeatedly differentiating the differential equation (j6.50p for its 
probability generating function and then letting z 1. In the supercritical 
case, this yields the moments of |Xak|1[|Xak| < oo] (or, equivalently, the 
moments of \Xxk\ conditioned on |Xak| < oo); it follows that all these mo- 
ments are finite, and we can obtain explicit formulae for them one by one. 
For example, with p = p{\k), 

E(|Xa.|2; < oo) = G"{1) + G'{1) = + i = i-, (6.52) 
E(|XA.h \^x.\ < oo) = G"'{1) + 3G"(1) + G'(l) = ^ + j-. (6.53) 

It can be seen that for each m > 1, as A \ Acr, and thus /? — > 0, we have 

E(|XA«r; I^AkI < oo) ~ cmp'-"" (6.54) 

for some constant Cm > 0; we do not know any general formula for Cm- For 
any A > Acr = | and a,b > 0, writing X := |Xak|1[|^Ak| < oo], from (j6.5ip 
and ([632]) - ([63l]) we obtain 

E(x2;X<^)<^EX = ^£(AK) = f, 
V pJ p p Xp 

k(x'^;X >-) < ^EX^ = + ^ 

and hence 



V ' - pJ - b hX^p bX' 



E(x^;^<X<^-)>^-^-^-^ = ^(l-a-^-P 
\ p pJ Xp Xp bX'^p bX Xp\ bX b 

Choosing, for example, a = 1/4 and b = 32, so 6A > 8, the last quantity is 
at least 1/(3A/?) > l.3/p if A is close to Acr, and thus, for such A at least. 

Hence, |Xak| may be as large as about p^^ with probability about p, as 
suggested by (|6.54p . 

Note that each p^^Xk) is a continuous function of A, so as A \ Acr, the 
(defective) distribution of |Xak| converges to the distribution of the critical 
|Xa„k|) which has mean x(-^cr'^) = 2 and P(|Xa„k| = A;) ~ 2/(A;^ log fc) as 
k ^ oo, see [17, Section 7.3]. 

In the subcritical case, PkiXK.) decreases as a power of k, see [13, Section 
7.3] for details. 
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We have so far studied xl-^^) ^-nd or, equivalently, the cluster size 

in the branching process Xak- Let us now return to the random graphs; we 
then have to be careful with the precise definitions. The Poisson version 
of the CHKNS model mentioned above can be described as the random 
multigraph where the number of edges between vertices i and j is Po(Ajj) 
with intensity Ajj := A(l/(7' — 1) — for 1 < i < j < n, independently 

for all such pairs i,j, see [igI : ItI : 3]- For the moment, let us call this random 
graph G^. Let G^/ be defined similarly, but with Xij := A(l/j — and let 
G^/^ be defined similarly with X^j := A(l/j — l/(n + 1)), for 1 < i < j < n. 
Since multiple edges do not matter for the components, we may as well 
consider the corresponding simple graphs with multiple edges coalesced; then 
the probability of an edge between i and j, i < j, is pij := 1 — exp(— Aij). (If, 
for simplicity, we consider A < 1 only, it is easy to see that the results below 
hold also if we instead let the edges appear with probabilities pij = Xij] this 
follows by the same arguments or by contiguity and [2Q . Corollary 2.12(iii)].) 

We first consider G^^; note that this is exactly (the Poisson version of) 
G^(n, Xk) with k defined in (j6.26p and the vertex space V given by 5 = (0, 1] 
with /X Lebesgue measure as above, and the deterministic sequence x„ = 
(xi, . . . , Xn) with Xi = i/n. Arguing as in the proof of Theorem l4.7l summing 
over distinct indices only, and using the fact that k is non-increasing in each 
variable, we find that the expected number KPi{Gl/) of paths of length i is 

^pdGi')< j2 n^^%^ 

jo,---,i£=i «=i 

— ^ / AK(xi_i, Xj) dxo • • • dxf 

j0,...,ji = l JU^iU^~^)/n,j^/n] -^^ 



< n 



\ TT A«;(xi_i,Xi) dxo • • • dx^ = n(rf^l, 1). 



Hence Lemmas 14.61 and 14.31 implv that (j4.5p holds and xiGn^) — ^ x('^'^)- 

For G^/^ , we observe that G^^^ can be seen as an induced subgraph of 
GI^j^i, and thus 

Ej^P,(G^^O <EE^^(^n+i) < (n + l)x(AK). (6.55) 

Hence Lemma HTSl imphes that xiG^^) — ^ x(Ak). 

Finally, it is easily checked that and G(/^ satisfy the conditions of (20I . 
Corollary 2.12(iii)], and thus are contiguous. Hence x{Gn) — ^ x(A«^) too. 
(One can also compare G^ and G^^ as in [3, Lemma 11].) 

It turns out that in probability bounds such as the one we have just proved 
do not obviously transfer from G^ to the original CHKNS model. On the 
other hand (as we shall see below), bounds on the expected number of paths 
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do. Hence, in order to analyze the original CHKNS model, we shall need to 
show that 

limsupEn-^^P^(G^) < x(Ak). (6.56) 

If A > 1/4, then \k supercritical, so x(Ak) = cxd and there is nothing to 
prove. Suppose then that A < 1/4. We may regard with the vertex 1 
deleted as GUIi- Writing P{G) for the total number of paths in a graph G, 
and P* for the number involving the vertex 1, by (j6.55p we thus have 

EP(G^)-EP*(G^) =EP(G^{_^i) <nx(AK), 

so to prove (1636)1 it suffices to show that EP*(G^) = o{n). 

Let S{GI^) denote the number of paths in G^ starting at vertex 1. Since a 
path visiting vertex 1 may be viewed as the edge disjoint union of two paths 
starting there, and edges are present independently, we have EP*(G^) < 
(ES'(G^))^. Now E5(G^) is given by 1 plus the sum over i of 1/i times 
the expected number of paths in Gj^Ii starting at vertex i. Durrett [l6l . 
Theorem 6] proved the upper bound 

3 1 (logz + 2)(logn - logj + 2) 
8 y/ij log n + 4 

on the expected number of paths between vertices i and j in the graph H 
on [n] in which edges are present independently and the probability of an 
edge ij, i < j, is l/(4j) (a form of Dubin's model; see the next section). In 
fact, his result is stated for the probability that a path is present, but the 
proof bounds the expected number of paths. (The factor l/y/Tj is omitted 
in (l6l . Theorem 6]; this is simply a typographical error.) This bound carries 
over to Gl^li, which we may regard as a subgraph of H. Multiplying by 
1/i and summing, a little calculation shows that this bound implies that 
ES'(G^) = 0(n^/^/ log n) for A = 1/4, and hence for any A < 1/4. From the 
comments above, (|6.56p follows, and for any A > we have x{Gn) — ^ x(Ak). 

Recall that the original CHKNS model G„ has the same expected edge 
densities as G^, but the mode of addition is slightly different, with or 1 
edges added at each step, rather than a Poisson number; this introduces 
some dependence between edges. However, as noted in |3], the form of this 
dependence is such that conditioning on a certain set of edges being present 
can only reduce the probability that another given edge is present. Thus, 
any given path is at most as likely in G„ as in G^ , and (j6.56p carries over to 
the CHKNS model. On the other hand, the effect of this dependence is small 
except for the first few vertices, and it is easy to see that Nk{Gn) has almost 
the same distribution as Nk{G(^). In particular, Nk{Gn)/n pkiXn), so 
the proof of Theorem 14.11 goes though. Using Lemma 14.21 it follows that 
xiGn) x(Ak). 

Turning to the supercritical case, let Mk{G) denote the number of com- 
ponents of a graph G, other than Ci, that have order k. We claim that, in 
all variants G^, G^,^ , G^^ or the original CHKNS model, for fixed A > Acr 
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there is some sequence of events £n that holds whp, and some r] > Q such 
that 

n-iK(Mfc(G„) I Sn) < lOOe-''*^'^', (6.57) 
say, for all n, > 1. Suppose for the moment that (16.571) holds. Then 

k>l k 

For each fixed k we have n-^¥.k'^Mk{Gn) = n"^ E(A;iVfc(G„) - 0{k)) 
kpk{XK). Since £n holds whp and n~^k'^Mi^{Gn) is bounded it follows 
that n^-*^/;;^ E(Mfc(G„) | £n) kpi.[\K). Hence, by dominated conver- 
gence, E(x(G'„) I £n) J2 ^Pki^f^) = xC-^/^)) and (which we know al- 
ready in this case), xiXf^) is finite. By Lemma l4.4( ii). it then follows that 

x{Gn) x(Ak) . 

To prove (I6.57|) we use an idea from [3]; with an eye to the next subsection, 
in the proof we shall not rely on the exact values of the edge probabilities, 
only on certain bounds. Fix A > Acr. Choosing ij small, in proving ()6.57p 
we may and shall assume that k is at least some constant that may depend 
on A. Set 6 = k~^^^^^, and let G'^ be the subgraph of Gn induced by 
the first n' = {1 — 6)n vertices. (We ignore the irrelevant rounding to 
integers.) In all variants G^, G^^, G^^^^ , the distribution of stochastically 
dominates that of Gn', so whp G^ contains a component G of order at least 
3p{XK)n' /A > p{\K)n/2. Let us condition on G^, assuming that this holds. 
Note that whp the largest component of G„ will contain G, so it suffices 
to bound the expectation of M^, the number of A;- vertex components of G„ 
not containing G. To adapt what follows to the original CHKNS model, we 
should instead condition on the edges added by time n' as the graph grows; 
we omit the details. 

Suppose that C is a component of G^ other than G. Consider some 
vertex v, n' < v < (1 — 5/2)n. Then v has probability at least X{l/v — 
1/n) > \5/{2n) > 5/{8n) of sending an edge to any given vertex, and hence 
probability at least 5\S\/(9n) of sending at least one edge to any given set 
S of vertices. Hence with probability at least (5^p(AK)|G'|/(200n), v sends 
an edge to both G and C . Since these events are independent for different 
V, the probability that C is not part of the 83/1116 component of Gn as G is 
at most 

(1 - (5V(AK)|G'|/(200n))'^"/^ < exp(-,5=^/5(AK)|G'|/400) = exp{-a5^\C'\), 

for some a > independent of k. 

Let A be the number of components of G'„ of size at least k^^^ that are 
not joined to G in G„. Then it follows that Ej4 < ne~"'^^ ^ . 

For any v < n' , the expected number of edges from 'late' vertices w > n' to 
V is at most 1/2, say. (We may assume 6 is small if A is large.) Let B be the 
number of vertices receiving at least /c^/^ edges from late vertices. Then it is 
easy to check (using a Chernoff bound or directly) that KB < ne^^^^^* for 
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some b > 0. The subgraph of G„ induced by the late vertices is dominated 
by an Erdos-Renyi random graph with average degree at most 1/2. Let 
N be the number of components of this subgraph with size at least k^^^. 
Then, since the component exploration process is dominated by a subcritical 
branching process, we have EA^ < ne~'^^^^* for some c > 0. 

Let M'f! be the number of fc-vertex components of Gn other than that 
containing C that do not contain any of the components/vertices counted 
by A, B or N. Since E(M^-M^') < E(A+S+iV) < ne-*^^'^" for some d > 0, 
it suffices to bound EM^'. Condition on G[j and explore from some vertex 
not in C. To uncover a component counted by M^', this exploration must 
cross from late to early vertices at least k^/^ times - each time we reach a 
component of size at most fc^^^, and from each of these vertices we get back 
to at most k^/^ late vertices, and from each of those to at most k^^^ other 
late vertices before we next cross over to early vertices. However, every time 
we find an edge from a late to an early vertex (conditioning on the presence 
of such an edge but not its destination early vertex), we have probability at 
least p{\k)/2 of hitting C. It follows that EM^' < n(l - p{\K)/2f^'\ and 
(16371) follows. 

Note that since x(^^) is a discontinuous function of A, we cannot obtain 
convergence to x(-^^) foi' an arbitrary sequence A„ A, as in Theorem 14.71 
and Section [6. li In fact, it follows easily from Theorem 14.11 that if A„ \ Acr 
slowly enough, then x(G'^("'5 A„k)) oo > x{^ctk) and x(G^(n,A„K)) > 
limA\Acr x(-^^) — e = 4 — e > x(Acr«^) whp for every e e (0, 2), for any vertex 
space V (with S and as above), and thus in particular for Gl^ . 

6.4. Dubins' model. A random graph closely related to the CHKNS model 
is the graph G^(n, A/t) with kernel 

K(x,y):=4- (6-58) 
X V y 

on 5 = (0, 1] , where the vertex space V is as in Section 16.31 so x„ = 
(xi, . . . ^Xn)- In this case, the probability pij of an edge between i and j 
is given (for A < 1) by pij = \K{i/n,j /n)/n = V j). Note that this is 
independent of n, so we may regard G^(n, \k) as an induced subgraph of an 
infinite random graph with vertex set N and these edge probabilities, with 
independent edges. 

This infinite random graph was introduced by Dubins, who asked when 
it is a.s. connected. Shepp [sO] proved that this holds if and only if A > 
1/4. The finite random graph G^{n,\K) was studied by Durrett [H, [13], 



who showed that Acr = 1/4; thus the critical value for the emergence of a 
giant component in the finite version coincides with the critical value for 
connectedness of the infinite version. See also (3: [29I: I3| . 
We have 

Tx.f{x) = - r f{y) dy + X I'M. dy. (6.59) 

X Jo Jx y 
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We can solve (I3.5p as in Section 16.31 we get the same equation ()6.28p 
and thus the same solutions f{x) = C4.x"+~^ + C_x"~~^ (unless A = 1/4 
when we also get a logarithmic term), and substitution into (I6.59P shows 
that this is a solution of (|3.5p if and only if C_|_q;+ + C_Q!_ = 1, see (j6.62p 
below. If < A < 1/4, so a+ > 1/2 > a_, there is thus a positive solution 
f{x) = a^^x°+~^ in L^. (This is the unique solution in L^, by a direct check 
or by Corollary 13.61 ) Hence, Corollary 13.61 yields 

x(Ak) = / fix) dx = = ^"^^ . < A < 1/4. (6.60) 



2A2 

Since this function is analytic on (0,1/4) but has a singularity at A = 1/4 
(although it remains finite there). Theorem 15.11 shows that Acr = 1/4, which 
gives a new proof of this result by Durrett [l^. Note that x('^cr'^) = 4 is 
finite. 

We can estimate the expected number of paths as in Section 16. 3^ and 
show by Lemmas 14.61 and 14.31 that x(G^('^7 -^^)) — ^ x('^'^) for A > 0. 
In the supercritical case, the tail bound ()6.57p goes through, showing 

that for any A > Acr we have xi^f^) < co, and x(G'^(™; -^^)) x('^'^)- 
Unfortunately, while the argument gives a tail bound on the sum kpi^(XK) 
for each fixed A > Acr, the dependence on A is rather bad, so it does not 
seem to tell us anything about the behaviour of xi^f^) ^ ^ approaches the 
critical point. 

We can easily calculate pk for small k. First, by ()6.59p . Txk^{x) = A — 
Alogx. Hence ()6.ip yields 

Pi{Xk; x) = e-^+^i°s^ = e'^x^. (6.61) 

Further, instead of ()6.3ip we now have, for every non-zero 7 > —1, 

Ta«(x^) = - - ^^x^ (6.62) 
7 7(7 + 1) 



Hence ()6.2p yields 

P2(Ak; x) = e~^x^Tx^{e-^x^) = e'^^x^ ( 1 - ) . (6.63) 

Similarly, by (|6.3p and some calculations, 

P-3A 

P3{Xk;x) 



x^ 



A + 1 



2(1 + A)2(1 + 2A) 
(2 + 3X)x^^ - 4(1 + 2A)(1 + X)x^^ + (2 + 3A)(1 + 2A)(1 + A)x^ 
and so on. By integration we then obtain 

/^i^^^) = TVx' ^^-^^^ 

^ 2Ae-2^ , , 

^^(^"^ = (1 + A)(1 + 2A)' ^'-''^ 
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, (15A2 + 18A3)e^3^ , , 

^^(^"^ = 2(1 + A)2(1 + 2A)(1 + 3A)- 

It is clear that each ^^(Ak) is e~^^ times a rational funtion of A, but we 
do not know any general formula or a recursion that enables us to calculate 
x(Ak) in the supercritical case as in Section lOl 

6.5. Functions of maxjx,?/}. The examples in Sections 16.31 and 16.41 are 

both of the type y) = ^pixN y) for some function Lp on (0, 1]. It is known 
that if, for example, '^(x) = 0(l/x), then is bounded on L^, and thus 
there exists a positive Acr > 0; see |2i;[l[ and 0, Section 16.6]. 
We have 

T^J{x) = \^{x) r f{y) dy + X [' ip{y)f{y) dy. (6.67) 



Jo Jx 

If (/? G C^(0, 1], then any integrable solution of (13. Sp must be in C^(0, 1] 
too, and differentiation yields /' = X(p'F, where F{x) := Jq f{y) dy is the 
primitive function of /; furthermore, we have /(I) = 1 + Tak/(1) = 1 + 
\ip{l)F{l). Hence, solving (jS.Sh is equivalent to solving the Sturm-Liouville 
problem 

F"{x) = Xip'{x)F{x) (6.68) 

with the boundary conditions 

F(0) = and F'(l) = Av9(l)F(l) + 1. (6.69) 

If there is a solution to (^6M\i and (fOQ]) with F' > and F' e L^, then 
Corollary 13.61 shows that 

x{Xk)= [ F'{x)dx = F{1). (6.70) 



JQ 

The examples in Sections 16. 31 and 16. 4l are examples of this, as is the Erdos- 
Renyi case in Section l6.ll (93 = 1). We consider one more simple explicit 
example. 

Example 6.7. Let (p{x) = 1 — x. Then ()6.68p becomes F" = —XF, with 
the solution, using (inSSD, F{x) = Asm{VXx) with A^/Xcos{^/X) = 1. This 
solution satisfies F' > if ^/A < tt/2, so we find Acr = and, by (j6.70p . 

x(Ak) = ^^^, A<A„ = vrV4. (6.71) 

6.6. Further examples. We give also a couple of counterexamples. 

Example 6.8. Let S = {1,2}, with fx{l} = fi{2} = 1/2, and let ^^(1, 1) = 
2, Ke{2, 2) = 1 and ^^(1, 2) = Ks{2, 1) = e for e > 0. 

For e = 0, kq is reducible; given the numbers rii and n2 of vertices of the 
two types, the random graph G^{n, Xkq) consists of two disjoint independent 
random graphs G{ni,2X/n) and G(n2,A/n); since ni/n,n2/n 1/2, the 
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first part has a threshold at A = 1 and the second a threshold at A = 2. Sim- 
ilarly, the branching process Xakq i^) is a single-type Galton- Watson process 
with offspring distribution Po(A) if x = 1 and Po(A/2) if j; = 2, so Xakq is 
a mixture of these. Hence, if XiW denotes the (modified) susceptibility in 
the Erdos-Renyi case, given by (|6.5p for A < 1 and (16. 6p for A > 1, then 

x(AKo) = ixi(A) + ixi(A/2), (6.72) 

so x(A/^o) has two singularities, at A = 1 and A = 2. Clearly, Acr = 1. 

Now consider e > and let e \ 0. Then Xcrii^e) < Acr(Ko) = 1- Fur- 
thermore, for any fixed A, p{Xks,x) —>■ p{\kq,x) by [1, Theorem 6.4(ii)], 
and hence — > (we may regard the operators as 2 x 2 matri- 

ces). Consequently, if A > 1 with A 7^ 2 and thus < 1; then 

(/ — T^)~^ {I — T^)^-*^, and thus xiX^^e) x(Ako) by Theorem 13.31 
This holds for A = 2 also, with the limit £(2^0 ) = ^1 for example by (|3.4|) 
and Fatou's lemma. 

Since x(A«^o) has singularities both at 1 and 2, we may choose 6 £ (0, 1/2) 
such that x((l + 5)'^o) > x(i«^o) and x((2 — (5)ko) > x(i'^o)) and then choose 
e > such that x((l + ^)i^e) > xi^i^e) and x((2 - S)ks) > xi^i^e)- This 
yields an example of an irreducible kernel k such that x(Ak) is not monotone 
decreasing on (Acr,cxD). 

Example 6.9. Theorem 14.71 shows convergence of x(G^(n, k)) to xi^^) for 
any vertex space V when k is bounded. For unbounded k, some restriction on 
the vertex space is necessary. (Cf. Theorem 14. 81 with a very strong condition 
on V and none on k.) The reason is that our conditions on V are weak and 
do not notice sets of vertices of order o(n), but such sets can mess up x- 

In fact, assume that k is unbounded. For each n > 16, find (a„, bn) G 5^ 
with K{an,bn) > n. Define x„ by taking [n^/^J points Xi = an, 
points Xi = bn, and the remaining n — points i.i.d. at random with 

distribution fi. It is easily seen that this yields a vertex space V, and that 
we have created a component with at least 2[n'^/^J vertices. Consequently, 
|Ci| > n3/4, and by x{G^{n,K)) > |Ci|Vn > n^/^, so x{G^{n,K)) 

00, even if k is sub critical and thus x('<^) < 

Using a similar construction (but this time for more specific kernels k), 
it is easy to give examples of unbounded supercritical k with x('^) < 00 but 
x{G^{n, k)) — > 00 for suitable vertex spaces V. 
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